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1.  Introduction 

Season  and  latitude  have  traditionally  featured  as  major  dependences 
of  the  structure  of  the  stratosphere  and  mesosphere.  The  Air  Force 
Reference  Atmospheres  (Cole  and  Kantor,  1978)  tabulated  temperature, 
density  and  pressure  for  each  month  and  each  15°  N  latitude  to  90  km 
altitude.  CIRA  1972,  Part  2  tabulated  these  parameters  for  each  month 
and  each  10°  N  latitude  to  110  km.  Tabulations  to  80  km  at  monthly 
intervals  have  been  extended  to  the  S  hemisphere  with  a  10°  latitude 
interval  (Koshelkov,  1983;  Barnett,  1984). 

Additional  dependences,  notably  local  solar  time  and  solar  activity, 
are  essential  for  modelling  of  the  thermosphere  as,  for  example,  in 
CIRA  1972,  Part  3  which  extends  upwards  from  110  km.  Further  dependences 
are  included  in  MSIS-83  (Hedin,  1983)  which  extends  upwards  from  90  km. 

The  present  report  addresses  itself  to  the  problem  of  modelling  the 
intermediate  region  between  mesospheric  and  thermospheric  models.  It 
seeks  to  develop  a  procedure  for  generating  a  smooth  transition  between 
selected  lower  and  upper  models  in  all  relevant  neutral  atmosphere 
parameters  including  composition. 

The  need  for  improved  modelling  of  the  mesosphere/thermosphere 
region  has  long  been  evident  and  a  1984  recommendation  for  a  new  CIRA 
referred  to  the  need  for  tabulations  from  80  to  120  km.  An  altitude 
range  from  70  to  130  km  has  been  adopted  for  the  present  work  to  allow 
a  smooth  transition  to  be  introduced  to  the  lower  and  upper  models, 
whose  values  are  required  to  be  matched  at  70  and  130  km.  The  problem 
of  modelling  the  70-  to  130  km  region  is  therefore  one  of  satisfying  two 
types  of  conditions.  One  is  imposed  by  the  requirements  of  the 
upper  and  lower  models  and  the  other  by  the  temperature  data  of  the 
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intermediate  region.  Such  data  are  taken  in  the  present  work  at  5  km 
height  intervals,  i.e.  at  75,  80,..., 125  km.  The  problem  is  a  novel  one 
in  atmospheric  modelling  in  that  the  usual  requirement  is  of  the  latter 
type,  i.e.  to  obtain  a  fit  to  a  set  of  data.  In  the  present  problem, 
the  requirement  to  also  match  lower  and  upper  models  in  all  their 
relevant  parameters  is  not  an  independent  one  and  a  particular  matter  of 
interest  and  concern  is  the  extent  to  which  both  types  of  condition  can 
be  satisfactorily  met.  The  treatment  here  is  to  rigidly  match  the  lower 
and  upper  models  (with  continuity  in  the  second  height  derivative)  and  at 
the  same  time  optimise  the  fit  to  the  intermediate  temperature  data, 
which  may  show  a  residual  bias  due  to  the  rigid  end  constraints.  In  prev¬ 
ious  work  (Forbes  and  Groves,  1986),  data  were  found  to  be  biassed  higher 
than  model  values  and  the  matter  to  warrant  further  consideration. 

The  procedure  of  model  formulation  is  described  in  Appendix  A  and 
relates  to  the  choice  of  MSIS-86  (Hedin,  1986)  for  the  upper  model  and 
to  that  described  in  'A  Global  Reference  Atmosphere  From  18  to  80  km' 
(Groves,  1985)  for  the  lower  model.  A  report  of  the  earlier 
calculations  undertaken  at  a  preliminary  stage  of  this  work  with 
MSIS-83  as  the  upper  model  has  been  given  previously  (Forbes  and  Groves, 
1986).  The  subsequent  stages  then  proposed  for  the  formulation 
have  now  been  completed  and  are  reported  here.  These  include 
(i)  modelling  of  individual  constituent  gases  and  hence  of  total  density 
and  pressure,  (ii)  combining  the  three  models  for  the  height  regions 
18-70,  70-130  and  above  130  km  to  provide  a  single  model  from  18  km  upwards 
(Appendix  D,  Section  3),  and  (iii)  the  development  of  models  for 
various  solar  conditions  (and  not  only  mean  solar  conditions  as  treated 
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earlier) . 


2.  Temperature  data  utilised,  70  to  130  km 

Temperature  data  have  been  assembled  from  two  earlier  collations  of 
data,  i.e.  from  that  used  for  the  construction  of  CIRA  1972,  Part  2 
and  the  rocket  and  incoherent  scatter  (I.S.)  data  reviewed  by  Forbes 
(1984).  Analyses  of  I.S.  data  from  St.  Santin  (Alcaycte  et  al.,  1979) 
have  also  been  utilised. 

I.S.  temperatures  were  utilised  as  monthly  means  for  the  locations 
and  heights  shown  in  Table  1  and  were  compared  with  models  computed  for 
the  solar  activity  parameters  shown  in  the  table. 

Table  1.  Distribution  of  available  monthly  mean  I.S.  temperatures. 


Site 

Height 

km 

Solar  activity 
F10.7=F10.7  Ap 

Millstone  Hill 

(42  N)  (Wand  1983,  Table  1) 

105-125 

120 

7 

St.  Santin  (45  N)  (AlcaydA  et  al.,  1979, 

Table  1) 

95-110, 120 

70,  120,  170 

7 

Arecibo  (18  N) 

(1970-75  data,  Forbes  1984) 

105-130 

108 

10 

Monthly  mean  temperatures  from  rocket  techniques  were  utilised  for  the 
locations  and  heights  shown  in  Table  2.  The  Kwajalein  data  of  1976-78  have 
been  compared  with  models  computed  for  F^q  -j  =  F^g  ^  =  95 units.  For  the 
other  data,  which  extend  no  higher  than  100  km,  solar  effects  are  minimal 
and  a  value  of  Fjq  y=  F^q  -j  =120  units  has  been  adopted  for  computing  models 
for  comparison  purposes.  Likewise  a  mean  value  of  (=10)  has  been 

adopted  in  model  computations . 


Table  2.  Distribution  of  monthly  mean  temperatures  obtained 
by  rocket  techniques. 


Site 

Height 

km 

Site 

Height 

km 

Heiss  Is.  (81  N) 

75-100 

Thumba  (8  N) 

75 

Volgograd  (48  N) 

75-100 

Ships  (equator) 

75-100 

Kwajalein  (9  N) 

75-120 

Molodezhnaya  (68  N) 

75-80 

The  locations  at  which  single  rocket  profiles  of  temperature  were 
utilised  are  shown  in  Table  3.  The  data  available  fall  off  rapidly  above 
100  km  as  shown  in  Table  3  by  the  numbers  of  profiles  available  at  105  and 
120  km.  Single  profiles  are  compared  with  models  computed  for  the  same 
solar  activity  parameters  as  defined  by  the  procedure  described  in  Appendix  B. 

Table  3.  Numbers  of  available  temperature  data  at  105  and  120  km 
from  rocket  launchings. 

Site  Height  (km)  Site  Height  (km)  Site  Height(km) 

105  120  105  120  105  120 

Pt  Barrow  (71 N)  3  0  Wallops  (38  N)  18  16  Barking  S  (22  N)  3  0 

Churchill  (59  N)  8  5  Whites  (32  N)  8  0  Kwajalein  (9  N)*  3  2 

Sardinia  (44  N)  0  1  Eglin  (30  N)  6  1  Ascension  (8  S)  2  0 

*  Developed  into  a  monthly  model  (Cole  et  al.,  1979)  and  listed  in  Table  2. 

3.  Model  fonrulation,  70  to  130  km 

The  procedure  devised  for  calculating  a  model  temperature  profile  from 
70  to  130  km  is  set  out,  step-by-step,  in  Appendix  A.  Theoretical 
details  are  presented  in  Appendices  Al  to  A5.  An  essential  feature 


of  the  method  is  the  expression,  for  given  geophysical  conditions,  of  g/T 
(g  being  gravity  acceleration  and  T  temperature)  as  a  polynomial  in 
height  whose  coefficients  are  chosen  (i)  to  give  continuity  up  to  the 
second  height  derivative  with  the  lower  and  upper  models,  (ii)  to 
reproduce,  on  integration  of  appropriate  physical  equations,  the  required 
ratio  of  ^  pressure  at  70  km  to  that  at  130  km  and  (iii)  to  produce  a 
best  fit  to  observed  temperatures  at  75,  80,..., 125  km  altitude. 

Appendix  A  sets  out  the  step-by-step  procedure  for  calculating 
density.  The  method  is  based  on  the  MSIS-86  (Hedin,  1986)  formulation 
for  number  densities  of  individual  gas  constituents.  Other  parameters 
then  follow  as  set  out  in  Appendix  A. 

4.  Comparison  of  observed  temperatures  with  models 

This  section  analyses  the  differences  between  observed  temperatures  and 
model  values  computed  in  each  case  for  the  same  geophysical  conditions. 
Model  values  are  derived  by  the  procedure  of  Appendix  A  which  involves 
polynomial  coefficients  agn  whose  determination  depends  on  the  analysis 
of  temperature  data  as  described  in  Appendix  B. 

Three  cases  with  different  selections  of  temperature  data  have  been 
considered  and  a  set  of  coefficients  agn  corresponding  to  the  second  of 
these  cases  is  listed  in  Appendix  B. 

4.1  Case  1:  a_  based  on  all  data  excluding  Millstone  Hill 
-  sn - - 

Case  1  was  undertaken  in  the  preliminary  investigation  (Forbes  and 
Groves,  1986)  which  utilised  all  data  at  first  and  then  excluded 
Millstone  Hill  data  because  of  their  relatively  high  values.  The 
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comparisons  reported  in  Sections  4.1.1  to  4.1.3  essentially  repeat  and 
update  the  results  previously  given. 

4.1.1  At  low  latitudes 

The  first  two  columns  of  Table  4  show  the  average  deviations  from  the 
computed  models  of  temperature  data  from  all  months  at  Kwajalein  (9  N) 
and  Arecibo  (18  N).  Such  deviations  are  denoted  by  x^  (for  the  i  th 
site)  in  Appendix  C  which  presents  the  relations  used  in  the  analysis. 
Kwajalein  and  Arecibo  provide  the  main  input  of  data  above  100  km  at  low 
latitude  and  as  previously  noted  (Forbes  and  Groves,  1986)  the  Kwajalein 
temperatures  are  (unexpectedly)  lower  than  those  of  Arecibo  by  about  15  K 
on  average. 

An  unsatisfactory  latitude  variation  at  low  latitudes  was  modelled 
in  the  preliminary  work  (Forbes  and  Groves,  1986)  by  closely  fitting  to 
these  data.  The  difficulty  has  been  overcome  in  the  present  work  by 
increasing  the  degree  of  smoothing  to  give  models  with  an  acceptable 
latitude  structure  at  low  latitude.  The  lower  Kwajalein  temperatures 
then  lead  to  lower  differences  as  shown  for  105  -  120  km  in  Table  4 r. 

In  the  last  column  of  Table  4-  are  the  mean  values  (x) ,  where  the 
mean  is  taken  over  all  sites  appropriately  weighted.  The  11  means  at 
75,...., 125  km  are  not  significantly  different  from  zero, 
having  an  average  of  1  K  (with  1  K  standard  deviation) .  AIK 
change  over  the  height  range  70  to  130  km  is  equivalent  to  a  4  per  cent 
change  in  the  ratio  of  ^  pressures  at  70  and  130  km.  Such  changes 
would  be  within  the  expected  limits  of  accuracy  and  therefore  the 
intermediate  temperature  data,  taken  as  a  whole  over  all  sites,  are  not 
inconsistent  with  the  lower  and  upper  models.  Data  for  a  particular 
site  may  nevertheless  still  be  biassed  one  way  or  the  other  with  respect 
to  the  models,  e.g.  Arecibo  data  at  105  to  120  km  show  a  region  of 
bias  to  higher  temperatures  by  about  10-  20 K. 


Table  4.  Low  latitude  temperature  differences  from  computed 
model  values  (K).  Av=  average  difference  for  data 
of  all  months  at  a  particular  site.  Case  1. 


Height  Kwajalein  Arecibo  Ships  All  sites* 

(equator) 


km 

Av 

sd 

Av 

sd 

Av 

sd 

Mean 
of  Av 

sd 

75 

-1.4 

1.3 

- 

- 

-0.5 

1.6 

-0.5 

1.7 

80 

2.4 

0.8 

- 

- 

-3.3 

1.2 

1.4 

1.5 

85 

5.9 

0.8 

- 

- 

-6.5 

1.0 

1.8 

3.2 

90 

6.9 

0.7 

- 

- 

-6.0 

1.2 

3.8 

2.8 

95 

3.9 

0.8 

- 

- 

-1.7 

2.1 

3.3 

1.7 

100 

-4.4 

0.9 

- 

- 

9.3 

4.6 

-3.5 

2.6 

105 

-8.0 

1.8 

9.9 

5.2 

- 

- 

-3.5 

13.0 

110 

1.4 

1.3 

13.9 

5.6 

- 

- 

2.1 

3.5 

115 

10.3 

3.9 

23.0 

6.1 

- 

- 

14.0 

8.2 

120 

-9.8 

2.7 

7.6 

4.6 

- 

- 

-5.3 

10.7 

125 

- 

- 

10.8 

7.0 

- 

- 

-10.8 

7.0** 

★ 

Ascension 

Is.  (8  S) 

,  Natal  (6  S) , 

Ships (equator) , 

Kourou  (5N) , 

Thumba(8N),  Kwajalein  (9  N) ,  Arecibo  ( 18  N) ,  Barking  Sands  (22  N) , 
Carnarvon  (25  S) . 
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Arecibo  data  only. 


4.1.2  At  middle  latitudes 

Table  5  presents  average  temperature  differences  from  computed  node Is  for 
the  5  middle  latitude  sites  that  contribute  most  data  at  heights  above  100  km.  • 
The  mean  of  such  averages  with  respect  to  all  sites  from  30°  to  50°  latitude 
is  also  shown.  In  contrast  to  the  results  of  Table  4,  model  values  are 
consistently  lower  than  observed  temperatures.  Such  a  bias  was  noted  and 
reported  at  the  preliminary  stage  of  this  work  (Forbes  and  Groves,  1986) 
and  is  able  to  arise  as  the  least-squares  fit  is  constrained  by  the  require¬ 
ment  to  match  the  ^  pressures  of  the  lower  and  upper  models  at  70  and  130  km. 

If  we  consider  the  changes  that  would  be  needed  in  the  models  at  70 
or  130  km  to  enable  the  intermediate  model  temperatures  to  fit  observed  values 
we  find  that  either  (i)  pressures  at  70  km  would  need  to  be  lower  by  40  per  cent 
or  (ii)  pressures  at  130  km  would  need  to  be  higher  by  40  per  cent  or  (iii) 
corresponding  partial  adjustments  would  need  to  be  made  at  both  heights. 

Such  pressure  adjustments  would  require  lower  temperatures  over  sane  range 
of  height  in  either  or  both  of  the  lower  and  upper  models  thereby  allowing 
higher  values  to  be  modelled  in  the  intermediate  region.  The  discrepancy 
between  observed  and  model  values. is  a  matter  of  discussion  in  Section  7  . 

4.1.3  At  high  latitudes 

Table  6  shows  average  temperature  differences  from  computed  models  for  three 
high  latitude  sites.  Data  are  available  from  few  sites  at  high  latitude 
and  then  mainly  below  105  km.  As  would  be  expected  in  these  circumstances, 
no  difficulty  arises  in  deriving  models  that  are  consistent  with  the  limited 
height  range  of  available  data. 

Table  6  shows  the  mean  of  the  average  differences  with  respect  to  the 
few  existing  high  latitude  sites.  The  absence  among  these  means  of  any 
significant  difference  from  zero  is  a  result  of  the  paucity  of  data  above 


Table  5.  Middle  latitude  temperature  differences  from  computed 
models  (K)..  Av  =  average  difference  for  data  of  all 
months  at  a  particular  site.  Case  1. 


Height 

Eglin 

White 

Sands 

Wallops 

St. 

Santin 

Mil 1st. 
Hill 

All  sites* 

km 

Av  sd 

Av  sd 

Av  sd 

Av  sd 

Av  sd 

Mean  of  Av  sd 

75 

-2 

2 

9 

4 

2  1 

. 

0.3 

1.1 

80 

1 

4 

9 

4 

0  2 

- 

- 

-  - 

1.1 

0.9 

85 

2 

3 

3 

5 

3  2 

- 

- 

.  _ 

5.8 

1.7 

90 

11 

3 

4 

4 

8  2 

- 

- 

-  - 

10.2 

1.9 

95 

20 

6 

15 

5 

14  4 

8 

1 

-  - 

9.8 

1.5 

100 

24 

7 

21 

7 

13  5 

1 

1 

-  - 

6.1 

5.1 

105 

24 

4 

27 

8 

5  5 

3 

1 

22  1 

13.2 

6.2 

110 

46 

19 

14 

8 

7  10 

10 

2 

38  3 

24.5 

12.0 

115 

- 

- 

- 

- 

54  14 

- 

- 

49  3 

48.9 

1.6 

120 

- 

- 

- 

- 

48  18 

32 

3 

39  3 

35.9 

4.6 

125 

- 

- 

- 

- 

27  25 

- 

- 

13  2 

12.8 

1.5 

* 


Eglin  (30  N) ,  Wocmera  (31  S) ,  White  Sands  (32  N),  Arenosillo  (37  N) , 
Wallops  (38  N) ,  Millstone  Hill  (42  N),  Sardinia  (44  N) ,  St.  Sant  in  (45  N) 
Volgograd  (48  N) ,  Kerguelen  (49  S). 


Table  6.  High  latitude  temperature  differences  from  computed 
models  (K).  Av  =  average  difference  for  data  of  all 
months  at  a  particular  site.  Case  1. 


Height 

km 

Churchill 

Av  sd 

Pt.  Barrow 

Av  sd 

Heiss 

Av 

Is. 

sd 

All  sites* 

Mean  of  Av  sd 

75 

3.6 

1.6 

-1.4 

2.2 

0.3 

1.3 

-1.4 

3.1 

80 

-0.2 

1.9 

-0.4 

2.4 

0.6 

1.8 

-1.1 

1.4 

85 

0.4 

1.7 

I 

o 

• 

3.0 

0.7 

1.6 

0.5 

0.4 

90 

5.8 

3.1 

6.7 

2.8 

-1.7 

2.0 

2.1 

4.9 

95 

11.2 

4.6 

15.7 

8.3 

-5.8 

2.7 

-0.2 

10.3 

100 

14.2 

4.1 

-7.6 

13.1 

-13.2 

2.8 

-4.6 

15.4 

105 

12.0 

7.4 

-24.8 

3.6 

- 

- 

-17.8 

20.5 

110 

6.0 

13.3 

- 

- 

- 

- 

6.0 

13.3** 

115 

26.3 

19.7 

- 

- 

- 

- 

26.3 

19.7** 

120 

47.3 

32.3 

- 

- 

- 

- 

47.3 

32.3** 

125 

- 

- 

- 

- 

- 

- 

- 

- 

*  Churchill  (59  N) ,  Molodezhnaya  (68  S) ,  Pt.  Barrow(7lN),  Heiss  Is.  (81N). 


**  Churchill  data  only. 


4.2  Case  2:  a  _  determined  without  I.S.  data  and  without  data 

- sn - 

from  105  -  125  km 

By  excluding  data  above  100  km  from  the  analysis  of  Appendix  B  and  the  determ¬ 
ination  of  agn,  we  can  expect  to  compute  models  that  have  a  better  fit  to 
data  at  and  below  100  km  at  the  expense  of  having  a  poorer  fit  to  the 
excluded  data  above  100  km.  Comparisons  between  observed  and  computed 
temperatures  are  presented  for  the  same  groups  of  low,  middle  and  high 
latitude  sites  as  for  Case  1.  Data  are  much  more  numerous  below  100  km 
than  above  and  for  this  reason  Case  2  is  worthy  of  investigation. 

4.2.1  At  low  latitudes 

Table  7  presents  the  results  for  Case  2  corresponding  to  those  for  Case  1 
in  Table  4.  As  would  be  expected,  model  temperatures  at  110-  125  km  are  now 
lower  (by  -  10  K)  and  those  at  85  -  100  km  are  higher  (by  ~  2  K).  These  changes 
are  generally  significant  for  any  particular  site  but  over  all  sites  the 
last  columns  of  Tables  4  and  7  show  that  the  changes  in  the  means  (which  are 
also  ~10K  and  2  K  in  the  respective  height  ranges)  are  not  significant, 
being  of  the  same  order  of  magnitude  as  the  sd's  of  the  means. 

The  fact  that  omission  of  data  above  100  km  results  in  such  a  limited 
change,  on  average,  in  the  modelling  accords  with  the  conclusion  of 
Section  4.1.1  that  observed  temperatures  for  75-  125  km  are  not  inconsistent 
with  the  lower  and  upper  models  when  taken  as  a  whole  over  all  sites. 

For  individual  sites,  the  average  differences  between  observed  and 
computed  temperatures  may  be  significantly  increased  or  decreased  for 
Case  2  relative  to  Case  1.  For  Arecibo,  observed  values  at  110-  120  km 
are  now  higher  than  computed  values  by  30-  40  K. 


-  11  - 


Table  7.  Low  latitude  temperature  differences  from  computed 
model  values  (K).  Av  =  average  difference  for  data 
of  all  months  at  a  particular  site.  Case  2. 


Height 

km 

Kwajalein 

Av  sd 

Arecibo 

Av  sd 

Ships 

(equator) 

Av  sd 

All  sites* 

Mean  sd 
of  Av 

75 

-1.3 

1.3 

- 

- 

-0.5 

1.6 

-0.4 

1.7 

80 

2.5 

0.8 

- 

- 

-3.7 

1.1 

1.1 

1.6 

85 

4.8 

0.7 

- 

- 

-7.5 

0.9 

0.4 

3.2 

90 

4.2 

0.5 

- 

- 

-7.8 

1.0 

1.7 

2.5 

95 

0.4 

0.8 

- 

- 

-3.6 

1.8 

-0.1 

1.4 

100 

-7.1 

1.0 

- 

- 

8.1 

4.3 

-6.2 

2.8 

105 

-8.0 

1.9 

9.3 

5.7 

- 

- 

-2.8 

14.9 

110 

5.4 

1.3 

23.3 

5.8 

- 

- 

6.4 

4.9 

115 

18.6 

4.0 

42.3 

6.2 

- 

- 

25.6 

15.3 

120 

-0.7 

2.7 

29.9 

5.3 

- 

- 

5.6 

17.5 

125 

- 

- 

-2.1 

6.8 

- 

- 

-2.1 

6.8** 

*  Ascension  Is. 

(8  S) , 

Natal 

(6  S), 

Ships  (equator). 

Kourou 

(5N) , 

Thumba  (8N),  Kwajalein  (9N),  Arecibo  (18 N),  Barking  Sands  (22  N), 
Carnarvon  (25  S). 

**  Arecibo  data  only. 


-  12  - 


4.2.2  At  middle  latitudes 

Table  8  presents  average  differences  for  Case  2  corresponding  to  Table  5 
for  Case  1.  As  expected  the  fit  up  to  100  km  is  much  improved  with  both 
positive  and  negative  average  differences  whose  means  over  all  sites  are 
now  not  significantly  different  from  zero  (the  last  column  of  Table  8). 

Above  100  km  the  high  average  differences  of  Table  5  become  still 
higher  for  Case  2  being  ^80K  at  115  -120  km. 

At  middle  latitudes.  Cases  1  and  2  present  straight  choices  between 
models  that  are  biassed  lower  than  observations  at  all  heights,  75-  125  km, 
(Case  1)  and  those  that  are  biassed  lower  than  observations  at  only 
105-  125  km  (Case  2),  the  biasses  nevertheless  being  significantly  greater. 

4.2.3  At  high  latitudes 

Table  9  shows  average  differences  at  high  latitude  sites  when  data  at 

105  km  and  above  are  omitted  from  the  determination  of  a  .  At  110  km 

sn 

and  above  the  omitted  data  amount  to  only  8  rocket  temperature  profiles 
at  Churchill  (59N)  and  hence  Table  9  shows  no  significant  change  from 
Table  6,  being  given  here  for  the  sake  of  completeness. 

4.3  Case  3:  a  _  determined  without  data  at  75  -  95  km 

- sn - 

In  this  case,  all  data  at  heights  100-  125  km,  including  incoherent  scatter 
data  from  Arecibo,  Millstone  Hill  and  St.  Santin  has  been  utilized  for 
determining  a  according  to  the  procedure  of  Appendix  B,  while  that  below 

Sli 
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Table  8.  Middle  latitude  temperature  differences  from  computed 
models  (K).  Av=  average  difference  for  data  of  all 
months  at  a  particular  site.  Case  2. 


Height 

Eglin 

White 

Sands 

Wallops 

St. 

Santin 

Millst. 

Hill 

All  sites* 

km 

Av  sd 

Av  sd 

Av  sd 

Av  sd 

Av  sd 

Mean  of  Av  sd 

75 

-1 

2 

10 

4 

3  1 

. 

1.0 

1.0 

80 

2 

4 

9 

4 

1  1 

- 

- 

_  _ 

2.0 

0.9 

85 

-2 

3 

-1 

5 

0  2 

- 

- 

-  - 

2.1 

1.8 

90 

-1 

3 

-6 

4 

-3  2 

- 

- 

-  - 

-0.3 

1.6 

95 

3 

5 

2 

4 

-1  4 

-6 

1 

-  - 

-2.9 

2.0 

100 

11 

7 

11 

6 

1  5 

-11 

1 

-  - 

-5.7 

5.3 

105 

23 

4 

26 

8 

4  5 

1 

1 

21  1 

13.8 

6.0 

110 

62 

19 

28 

9 

25  10 

23 

2 

53  2 

42.4 

12.8 

115 

- 

- 

- 

- 

91  14 

- 

- 

80  3 

80.8 

2.9 

120 

- 

- 

- 

- 

90  17 

65 

3 

75  3 

69.9 

6.3 

125 

- 

- 

- 

- 

45  25 

- 

- 

27  2 

27.2 

1.7 

* 


Eglin  (30N),  Wocmera  (31  S),  White  Sands  (32N),  Arenosillo  (37  N), 
Wallops  (38N),  Millstone  Hill  (42  N),  Sardinia  (44N),  St.  Santin  (45  N) 
Volgograd  (48 N),  Kerguelen  (49  S). 


Table  9  High  latitude  temperature  differences  from  computed 
models  (K).  Av  =  average  difference  for  data  of  all 
months  at  a  particular  site.  Case  2. 


Height 

km 

Churchill 

Av  sd 

Pt.  Barrow 

Av  sd 

Heiss 

Av 

Is. 

sd 

All  sites* 

Mean  of  Av  sd 

75 

4.3 

1.6 

-1.2 

2.2 

0.2 

1.3 

-1.2 

3.2 

80 

1.5 

1.9 

0.3 

2.4 

0.4 

1.8 

-0.3 

1.4 

85 

0.4 

1.7 

0.0 

3.0 

0.6 

1.6 

0.4 

0.2 

90 

1.9 

3.0 

5.9 

2.8 

-1.6 

2.0 

1.2 

3.9 

95 

3.9 

4.7 

14.2 

8.2 

-5.3 

2.7 

-1.7 

7.3 

100 

7.5 

4.2 

-9.3 

13.1 

-12.7 

2.8 

-6.6 

11.3 

105 

9.4 

7.6 

-25.7 

3.6 

- 

- 

-19.3 

19.2 

110 

10.8 

13.6 

- 

- 

- 

- 

10.8 

13.6** 

115 

40.0 

19.9 

- 

- 

- 

- 

40.0 

19.9** 

120 

64.3 

32.1 

- 

- 

- 

- 

64.3 

32.1** 

125 

— 

_ 

_ 

_ 

_ 

_ 

_ 

*  Churchill  (59  N) ,  Molodezhnaya  (68  S) ,  Pt.  Barrow(71N),  Heiss  Is.(8lN). 

**  Churchill  data  only. 
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100  Van  has  been  omitted.  In  comparison  with  Cases  1  and  2  an  improvement 
in  the  fit  of  computed  temperatures  to  observed  values  can  be  expected  above 
100  lan  at  the  expense  of  lower  computed  temperatures  and  a  poorer  fit  to 
observations  below  100  Van. 

Table  10  shows  the  results  obtained  for  mid-latitude  data  corresponding 
to  Tables  4  and  7  for  Cases  1  and  2.  At  85  -  95  Von,  Table  10  shows  that 
temperatures  are  '*'20  K  higher  than  model  values.  In  comparison,  the 
difference  is  '*'10 K  for  Case  1,  but  differences  above  100  km  are  greater. 
Case  2  produces  a  close  fit  at  85-  95  km  with  still  larger  differences  above 
100  km. 

5.  Tabulations  of  temperature,  pressure  and  density 
A  wide  choice  of  options  is  open  for  the  particular  models  to  present  as 
tables  and  for  the  format  in  which  they  are  to  appear.  The  choice 
involves  the  height  interval  of  tabulation  and  many  geophysical 
parameters . 

Diurnal  and  zonal  mean  values  have  been  tabulated  in  the  course 
of  this  project  as  latitudinal -height  cross-sections  for  solar  activity 
parameters  that  are  common  with  MSIS-86  tabulations,  i.e.  F^q  -j  =  70, 

150,230  units  and  A^  =  4,  48,  400. 


r 


E* 


I*. 


Table  10.  Middle  latitude  temperature  differences  from  computed 
models  (K).  Av=  average  difference  for  data  of  all 
months  at  a  particular  site.  Case  3. 


Height 

Eglin 

White 

Sands 

Wallops 

St. 

San tin 

Mil 1st. 
Hill 

All  sites* 

km 

Av  sd 

Av  sd 

Av  sd 

Av 

sd 

Av  sd 

Mean  of  Av  sd 

75 

0 

2 

12 

4 

5 

1 

_ 

2.6 

1.2 

80 

10 

4 

20 

4 

11 

2 

- 

- 

- 

- 

10.8 

1.0 

85 

16 

3 

19 

5 

19 

2 

- 

- 

- 

- 

20.3 

1.7 

90 

23 

3 

18 

4 

23 

2 

- 

- 

- 

- 

23.4 

1.9 

95 

25 

5 

21 

5 

22 

4 

16 

1 

- 

- 

18.5 

1.6 

100 

19 

7 

17 

7 

10 

5 

-1 

1 

- 

- 

4.4 

5.3 

105 

10 

4 

10 

8 

-12 

5 

-12 

1 

6 

2 

-6.5 

5.7 

110 

23 

18 

-14 

8 

-23 

9 

-19 

2 

8 

4 

-2.4 

14.1 

115 

- 

- 

- 

- 

19 

12 

- 

- 

13 

4 

13.3 

2.6 

120 

- 

- 

- 

- 

23 

16 

6 

4 

13 

4 

9.9 

5.3 

125 

- 

- 

- 

- 

22 

25 

- 

- 

6 

1 

5.9 

1.3 

*  Eglin  (30  N),  Woomera  (31  S) ,  White  Sands  (32  N),  Arenosillo  (37  N) , 

Wallops  (38  N) ,  Millstone  Hill  (42  N),  Sardinia  (44  N) ,  St.  Santin(45N), 
Volgograd  (48  N) ,  Kerguelen  (49  S) . 
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Two  formats  have  been  adopted  and  coded  (Appendix  D,  Section  3). 

One  has  a  5  km  height  interval  (70-  130  km)  and  20°  latitude  interval 

(80S  to  80 N)  with  12  such  tables  per  page  consisting  of  temperature, 

pressure  and  density  cross-sections  for  four  selected  months.  Tables 

are  presented  in  Appendix  F  for  all  months  based  on  Case  2  determinations  of 

a  and  solar  activity  F,n  =  150  units,  A  =4. 
sn  J  1U./  p 

The  second  format  has  a  1  km  height  interval  and  10°  latitude 
interval  from  80S  to  80 N  and  is  identical  with  the  format  previously 
used  (Groves,  1985)  for  18-80  km  with  each  table  occupying  one  page. 

The  height  range  is  from  65  to  135  km  with  the  values  at  65  -  70  km  being 
those  of  the  lower  model  and  at  130  - 135  km  being  those  of  the  upper  model 
i.e.  MSIS-86.  This  format  is  presented  in  Appendix  F  for  January  temperature, 
pressure  and  density  cross-sections  based  on  Case  2. 

6.  Limitations  of  the  model  formulation,  70-  130  km 

6.1  Types  of  data 

The  only  type  of  data  that  has  been  utilized  at  heights  75-  125  km  is 
temperature,  which  has  served  both  as  an  input  to  the  formulation  of  the 
models  for  these  heights  and  for  comparison  with  the  computed  models  to 
check  their  goodness -of -fit.  A  computed  composition  is  available  from  the 

fornulation,  but  compositional  data  have  not  been  utilized  either  as  an 
input  to  the  formulation  or  to  check  against  computed  compositions. 

6.2  Local  solar  time 

No  representation  of  tidal  components  at  70  -  130  km  is  included  in  the 
fornulation.  At  130  km,  dependence  on  L.S.T.  is  that  defined  by  MSIS-86 
while,  at  70  km,  the  lower  model  is  independent  of  L.S.T.  At  intermediate 
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heights  the  L.S.T.  dependence  arises  by  interpolation  vrtule  ignoring 
tidal  fields  that  may  have  detailed  spatial  structures. 

6.3  Longitude  and  solar  activity  effects 

No  detailed  representation  of  longitudinal  variations  at  70-  130  km  is 
included  in  the  formulation.  At  70  and  130  km,  the  dependences  on  longitude 
are  those  of  the  adjacent  lower  and  upper  models  and,  relative  to  other 
variations  at  these  heights,  are  very  small.  At  intermediate  heights, 
longitude  dependence  is  then  generated  by  interpolation  and  is  a  correspondingly 
small  effect. 

Dependence  on  solar  activity  is  also  not  directly  represented  at  heights 
between  70  and  130  km,  but  is  present  in  the  interpolation  between  these 
two  heights  where  it  is  represented. 

6.4  The  coefficients  a  _ 

- sn 

Sets  of  polynomial  coefficients  a  (s=l,...,S;  n  =  l,...,N)  are  introduced 
by  which  height- latitude  cross-sections  of  atmospheric  structure  parameters 
may  be  generated  from  70  to  130  km  at  all  latitudes.  Subscript  n  is  associated 
with  a  polynomial  in  height  and  subscript  s  with  a  polynomial  in  sin( latitude ) . 
Lack  of  data  limits  the  number  of  sets  of  coefficients  determined  to  12, 
one  for  each  month  of  the  year.  In  principle,  additional  sets  could  be 
introduced  with  dependences  on  L.S.T. ,  longitude  and  solar  activity  to  help 
to  overcome  the  limitations  mentioned  in  Sections  6.2  and  6.3,  but  a  vastly 
greater  amount  of  data  than  that  available  would  be  required. 

After  many  test  computations  with  different  choices  of  N  and  S,  the 
choice  was  made  of  N  =  2  and  S=7,  giving  14  coefficients  per  set.  Without 
improvements  in  the  quality  and  quantity  of  available  temperature  profiles  and 


their  latitude  distribution,  N  and  S  could  not  be  justifiably  increased. 

Since  agn  relate  to  a  pole-to-pole  representation,  they  provide  a 
means  of  formulating  seasonal  asymmetries  between  N  and  S  hemispheres. 
Unfortunately  observations  are  not  available  for  much  of  the  S  hemisphere 
and  it  becomes  necessary  to  assume  seasonal  symmetry  for  latitudes 
where  observations  are  lacking  in  the  determination  of  agn  (Appendix  B). 

7.  Discussion 

The  modelling  of  atmospheric  properties  in  an  intermediate  height  region 
(70-  130  km)  between  given  mesospheric  and  thermospheric  models  has  been 
treated  theoretically  (Appendices  A,  A1  to  A5  and  B)  and  computationally 
as  summarised  in  Appendices  D,  E  and  F. 

Temperature  data  (at  75,80,...,  125  km)  that  are  utilized  for  modelling 
the  intermediate  region  are  sumnarized  in  Section  2  and  two  formats,  that 
have  been  devised  for  tabulating  models  of  temperature,  pressure  and  density, 
are  outlined  in  Section  5  with  examples  given  in  Appendix  F. 

Section  4  presents  comparisons  between  observed  temperatures  (75-  125  km) 
and  model  values  and  is  a  Section  that  corrmands  particular  attention  on 
account  of  biasses  between  the  two  which  point  to  a  possible  inconsistency 
between  the  observed  temperatures  (75-  125  km)  and  the  given  lower  and  upper 
models  that  are  matched  at  70  and  130  km.  The  biasses  feature  strongly  at 
mid- latitudes,  where  the  majority  of  available  data  have  been  obtained,  and 
are  in  the  sense  that  observed  values  are  higher  than  model  values  (Case  1 
and  Table  5) . 

One  way  in  which  the  discrepancy  could  be  resolved  would  be  by  modifying 
the  N2  pressures  of  the  lower  and  upper  models  at  one  or  both  of  the  heights 
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70  and  130  km.  The  magnitude  of  the  discrepancy  is  such  that  its  removal 
would  require  a  40  per  cent  change  in  the  ^  pressure  ratio  for  these  two 
heights  and,  as  only  as  few  per  cent  adjustment  could  reasonably  be  allowed 
at  70  km,  the  majority  of  it  would  need  to  be  applied  at  130  km. 

Another  possible  conclusion  is  that  observed  temperatures  are  erroneously 
high  over  at  least  some  part  of  the  70-  130  km  height  region.  Above  100  km, 
temperature  data  are  provided  by  the  I.S.  technique  and  a  small  number  of 
measurements  by  rocket  techniques  as  summarised  in  Tables  1  to  3.  Up  to 
100  km  rocket  data  are  much  more  numerous  and  are  obtained  by  well-established 
techniques,  such  as  the  grenade  experiment.  Case  2  (Section  4.2)  was 
therefore  introduced  to  provide  models  that  are  consistent  with  the  data  at 
least  up  to  100  km  (Table  8).  The  closer  fit  up  to  100  km  is  however 
obtained  at  the  expense  of  a  poorer  fit  above  100  km  as  shown  in  Table  8, 
where  biasses  at  115-  120  km  are  close  to  80 K  for  both  I.S.  measurements 
at  Millstone  Hill  and  St.  Santin  as  well  as  rocket  techniques  at  Wallops. 

Support  for  high  I.S.  temperatures  is  provided  by  Arecibo  data  (18N) 
as  shown  in  Table  7,  where  biasses  at  115-  120  km  are  about  30-  40 K,  i.e. 
about  half  those  at  mid-latitudes .  The  presence  of  a  small  component  of 
high  energy  electrons  would  give  I.S.  temperatures  in  excess  of  neutral 
gas  temperatures  and  could  possibly  account  for  the  discrepancy.  The 
presence  of  a  similar  bias  in  the  rocket  measurements  at  Wallops  (Table  8), 
which  at  120  km  amount  to  16  in  number  (Table  3),  means  that  the  discrepency 
is  unlikely  to  be  solely  attributable  to  biassed  I.S.  temperatures.  In 
particular  the  possibility  of  adjusting  to  higher  ^  pressures  at  130  km 
should  be  considered  as  by  that  means  higher  model  temperatures  at  115-  120  km 
would  result  and  reduce  the  discrepancy  whatever  the  method  of  measurement. 


Finally,  the  contrary  view  to  that  represented  by  Case  2  might  be 
taken,  namely  that  the  inconsistency  is  attributable  to  observed 
temperatures  below  100  km  being  too  high.  Case  3  therefore  computes 
temperature  models  by  fitting  to  data  at  and  above  100  km  only  and  in  so 
doing  generates  values  at  85  -  95  km  that  are  on  average  ~20K  lower  than 
observed  values  (Table  10).  At  these  heights,  temperatures  have  been 
measured  by  the  grenade,  falling  sphere  and  pitot  pressure  techniques, 
which  are  well-established  methods  that  would  not  be  expected  to  be 
consistently  biassed  by  more  than  a  few  degrees  K  and  certainly  not 
by  as  much  '*-20  K.  Case  3  is  not  therefore  considered  to  provide 
an  acceptable  model. 

No  positive  recommendation  for  acceptable  models  of  the  intermediate 
(70-  130  km)  height  range  has  been  reached  in  the  above  investigations 
on  account  of  consistent  differences  between  temperature  data  and  the 
possible  models  that  can  be  computed  with  matching  conditions  at  70  and 
130  km.  One  source  of  such  differences  between  data  and  models  could 
be  the  method  of  model  formulation  and  the  possibility  of  incompletely 
formulated  dependences  on  L.S.T.,  longitude  and  solar  activity  is  pointed 
out  in  Section  6.  The  differences,  however,  appear  to  have  no  consistent 
relationship  with  these  dependences  and  are  much  greater  in  magnitude 
than  any  expected  shortcomings  in  the  formulation  with  respect  to  these 
dependences . 

Attention  therefore  needs  to  be  given  to  the  possibility  of  observed 
values  of  temperature  being  overestimated.  Case  2  has  been  devised  to 
fit  observed  temperatures  closely  up  to  100  km,  but  in  so  doing  it  provides 
models  above  100  km  that  are  less  than  observed  mid-latitude  values  by 


— ’  80 K  at  115  -120  km  (Table  8)  in  order  to  match  boundary  conditions  at 
70  and  130  km.  There  may  be  scope  for  reducing  the  ^80K  difference  at 
115-  120  km  by  revising  MSIS-86  at  130  km  -  the  region  from  130  to,  say, 

150  km  being  poorly  observed  -  but  the  general  nature  of  the  inconsistency, 
i.e.  observed  temperatures  in  excess  of  model  values,  would  still  remain. 
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APPENDIX  A 


MODEL  FORMULATION 


At  height  z  in  the  range  (ZpZ2),  where  z^  ^  80  km  (80  km  being  the  upper 
limit  of  the  lower  model  (1))  and  z2>90  km  (90  km  being  the  lower  limit 
of  the  upper  model  (2)),  we  write 


T-  (\t/R)  W 


y»-0 


where  W  is  a  function  of  height  whose  determination  is  the  central 


consideration  of  the  model  formulation.  g  is  acceleration  due  to  gravity, 

R  the  universal  gas  constant  and  M  the  mean  molecular  weight  of  air  as 

Z1 

given  by  the  lower  model  at  height  Zy 

As  R  may  have  (slightly)  different  values  for  any  given  lower  and 
upper  models  and,  likewise,  g  may  be  represented  by  different  expressions 
in  the  two  models,  the  values  of  R  and  g  at  height  z  are  based  on  a 
smooth  transition  from  R^,  gj(z)  for  the  lower  model  to  R2,  g2(z)  for  the 


upper  model  according  to  the  relations 


R  -  if  +  R,  -  cSl 

§  =  4  [  3,  *  U  +  CjrSOfc-K 


where 


0 


(A.-u) 


i 

(M) 


~8 

With  c  =  10,  the  accuracy  at  }  =  -1  is  better  than  1  in  10  .  The  lower 

and  upper  models  provide: 

M  *  HS.'jCk 4  (■&>»{> Jt)  R,~  8' 5*1*  2-?  *  10  ^  8.5 i4  >(o  JK 
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WoSS’Lil  +  O-OoSZMS’^Sf  -  o.ooooo^ ^x2f)  <>*"*)  (*•«) 


rf  -  /  (3.offr46  +  0-00T2jc^2<t>)  (JL.)  (a.<?) 

§z  =  *S  /  (  1  +  y/^T  CA-'°) 

where 

<JS=  1.30U5T  •»,**  f  635*6.77^  (Am) 

(It  may  be  noted  that,  if  the  mean  molecular  weight  of  air  is  constant  for 
z  •*£■  ,  W  is  pressure  scale  height  (km)  for  z  «  z^ ) . 

For  Zj  ss£  z  **■  ,  we  express 

vr'=  A  +  6  (^-‘)  (A.U') 

where  A  and  B  are  polynomials  in  5  .  B  is  an  interpolating  polynomial 
which  depends  only  on  conditions  at  height  z^  and  z^,  as  defined  by  the 
particular  lower  and  upper  models  under  consideration,  and  is  otherwise 
independent  of  conditions  relating  to  the  interval  .  The  conditions 

imposed  at  heights  z^  and  Z2  are  however  such  (see  below)  that  B  may  be 
considered  to  provide  a  first  approximation  to  W~*. 

A  is  an  ’adjusting*  polynomial  whose  determination  is  independent  of 
conditions  at  heights  z^  and  Z2  being  dependent  on  temperature  observations 
in  ( z^ , Z2 )  or  more  correctly  on  the  differences  between  values  of  W  ^ 
calculated  by  (A.l)  from  available  temperature  values  and  B  calculated  for 
the  same  geophysical  conditions  (of  date,  location,  solar  activity  etc.). 

We  take  A  and  its  first  two  height  derivatives  to  be  zero  at  heights 
Zj  and  Z2  by  expressing  it  as 

a-  (i-jfz  I  «.,h-r.)54'  ; 

*  I  n*  I 

We  choose  =  0,  =  1/9  as  explained  in  Appendix  Al.  The  determination 
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of  agn  is  described  in  Appendix  B,  where  values  of  agn  are  tabulated. 

B  is  determined  by  7  conditions  involving  model  values  at  heights  z^ 
and  z2  and  is  therefore  taken  as  a  polynomial  of  degree  6  in  5 

3  -  +  txS  -»• _ i7  *  [i  s  - 

where  b  =  [b^  b2  . . .  by .  Six  of  the  conditions  are  for  continuity  of 
W  (and  hence  of  T)  and  of  its  first  and  second  height 
derivatives  at  heights  z^  and  The  seventh  condition  is  that  the 

ratio  of  the  N2  pressures  at  and  z2  as  calculated  frcm  the  model 
temperature  profile  should  equal  that  of  the  N2  pressures  specified  by  the 
lower  and  upper  models  at  z^  and  z^.  We  write  (Appendix  Al) 
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For  the  conditions  at  height  z^ ,  we  have  (Appendix  A2) 

i,  -  4,  /  3>,  A,  (*•'») 
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A  =  (A/:D.’+  )>* 
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(A. 3a) 

-*r»  cr=i>j'3-) 

(A.Zi) 

Hre^  is  the  zonal  mean  pressure  scale  height  of  the  lower  model  and  hence 
from  (A. 24)  and  Ref.  1  we  have 


J~.AU/  H (Ms) 


-ft. 


n-l  yj»l 
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-fl  * '  s6wl 


(4-2)  (MO 


where 

Z, *  C%-2„)IZjl  ,  2.„=  48.75 -A,  Zj=V.lS-L  fA.2S) 

and  cng  are  linearly  interpolated  to  the  required  date  from  the  values 
tabulated  in  units  of  km  ^  in  Ref.  1.  The  dependence  of  ^.l^.l^  on 
longitude,  X  ,  is  expressed  by  q  ,  which  by  (A2.3)  and  (A 2.7)  is 


—  Klfm  Cr$  X  L)r  -4^  X  c‘r4  ^  X  l-^r  •<£*■>  2  X  C  A. 

vrtiere,  for  j  =  1,2;  r  =  1,2,3, 


-  29  - 


*jr  =  r.  'i 

*-jv 
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\  t~ 
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r  =  CA-5°^ 

»/ 


C&5 


If  ~  ^  Vy  /  10 

Values  of  y  are  calculated  from  the  values  of  Tj,  tabulated  in  Ref.  1 
for  each  month  at  10°  latitude  steps  and  are  linearly  interpolated  to  the 
required  date  and  latitude.  Yr  are  obtained  numerically  from  values  y_^, 

yi’y3/2’y5/2  °f  y  at  ^  equally- spaced  heights  z_|>z^>z3/2»z5/2  which  zi 
is  the  mid-point  and  Az  is  the  increment  step  by  fitting  a  cubic  to  provide 

Y,  =  &+)/«> 

Y.  =  (27a.— (^-53) 
Y3  =  («+-  t+) I  *  faf  (*-**>) 


where 


a±  =  ^  ~  ±  2-^ 

For  the  tabulations  in  Ref.  1,  fa z  =  4  km.  For  z^  =  70  km  (the  value 

adopted),  the  4  equally-spaced  heights  are  64,68,72  and  76  km.  (For  other 

values  of  z^,  (A. 32)  to  (A. 34)  may  need  to  be  replaced  by  alternative 

2/  j) 

formulae).  The  dimensions  of  qr  are  (km)  ^  for  r  =  1,2,3. 

The  terms  Ar  relate  to  incremental  changes  of  temperature  in  the 
vicinity  of  height  z^  that  may  be  associated  with  the  solar  cycle. 
Corresponding  to  a  change  ARn  in  sunspot  number  from  a  reference  value 


(A.Ss) 


where 


we  have  the  formulation  (Appendix  A2) 
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1,  =  ( !-*’)(/ +*«’*)  CA.») 

fz  =  tC,-t*X-2  +3*t  -  s-«t  «■*)/-£.  (a.1<j) 

where 

=  tarrJkltyro)lJt]  ^'Uy^ 

The  quantities  oC,a,l,p,q  and  n  are  calculated  from 

&  ~  4  6Z  (  9  =  ct,a,l,p,q  or  n)  (A.^a) 

^  —  tasoJh  (/4-^/lc)  c<^\yrC^'<L~'[)  I  (/^.Zi3) 

t^  being  the  day  number  in  the  year.  Numerical  values  for  6^ ,  ^ 

(  0  =  <^,a,l  ,p,q  or  n)  have  been  given  (3)  but  are  tentative  being  derived 
from  limited  data  at  eastern  longitudes  and  therefore  this  dependence  has 
only  been  included  in  exploratory  calculations  that  are  not  a  part  of  this  report. 
The  model  of  Ref.  1  is  for  data  over  years  of  average  sunspot  number,  RnQ=65. 

For  the  conditions  at  height  Z2,  we  have  (Appendix  A3) 

A  =  10,3  M$,  KTi*,)  (A-4^ 

where  T^)  is  the  temperature  at  height  Z2  calculated  from  MSIS-86  for  the 
required  time,  date,  latitude,  longitude,  solar  activity  etc.  Also 


Ca.^ 

-£7  =  (U  +  ^^+7)^3  ^ 

(a.m) 

where 

d  ~  CKj»  +  y2)  ‘ 

(-fc-m  ’  ) 

(A.47) 

For  z0  ^  z  =  117.2  km 

2  a 

04.  =  -  Sjj.  [  1  ~  /^To  )] 

(A. 48) 

•ft  *  1L  t2.lL  -*■  ) 

(A -4<^ 
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where  <T  is  taken  from  MSIS-86,  being  related  to  the  MSIS-86  parameters  of 
temperature  and  temperature  gradient  at  =  120  km  and  the  MSIS-86 
exospheric  temperature  by  S'  =  /  (T,*,  -T^ ) . 


For  z~^  z  =  117.2  km, 
L  a 


(A.S2) 


where 

**  =  -  r ? Or, lO  -  S’ £r„ ,i„ l]  /  ? ")  (*■ ? 

with 

5Ctt,x)  =  +  (*■* 

and 

/  4Zr\  -  —  1  |_^£)  (A.s 

The  coefficients  Tg,T^,T^  are  taken  from  MSIS-86  being  related  to  the 
MSIS-86  parameters  of  temperature  T  and  temperature  gradient  T’  at  z  =  z 
and  the  MSIS-86  temperature  Tq  at  height  zq  of  the  MSIS-86  mesopause. 
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Finally,  for  the  seventh  condition,  which  involves  values  at  both 


heights  Zj  and  z^  we  have  (Appendix  A4) 

A  =  (>•«) 

where  =  28  kg/(kmol)  ^  and  X  is  obtained  by  iteration  (over  only  2  or 
3  cycles  for  10  ^  accuracy  on  taking  an  initial  value  of  X  =  ^-z^/Okm)) 
from  (A4.9)  written  as 

x  -  CM*  /  ff.  )  lA  A-  <  ^  ( 1  +  A"'-X  ll  ^7) 

where  Mq  =  28.95  kg(kmol)-^  and 

^  (V,  ,  M,Vj  )  ^-5")) 

aKi=  ^I(mo-m^)  A**) 

where  =  28  kg(kmol)  .  is  the  fractional  volume  of  air  at 

height  z^  (=  0.78084)  and  p(z^),  the  air  pressure  at  height  z^,  is 

(ifyi)  =  (v,  )  C  I  -+'3?c)Cl  +Ao)  (A.6l) 

where  Pref(z^)  is  the  zonal  mean  air  pressure  at  height  z^  calculated 
from  Ref.  1 

1  1  *'*  V  \  . 

)WO()  “■  f-31.25'  I  i.  f  5/^  (^2) 

ip  /v-  n  =o 

( 5°/°  denoting  unity  for  n  =  0).  cng  are  linearly  interpolated  to  the 
required  date  from  the  values  tabulated  in  Ref.  1.  Dependence  of  p(z^ ) 
on  longitude  N  is  introduced  through 


K]Q  cd  \  "+■  l-to  ^  +  ^20  cn>  2  -t-  iL26  -du*  2.^  (A£Z^ 


where,  for  j  =  1,2, 


L  °  =  {^\ 

Uic  'I 


'  J=NZhi 


(*■£*) 


Values  of  y  are  calculated  from  the  values  of  p ^ ,  )\^  tabulated  in  Ref.  1 
for  each  month  and  10°  latitude  step  and  are  linearly  interpolated  to  the 
required  date  and  latitude.  The  value  of  y  at  height  is  obtained  by 
an  interpolating  cubic  through  four  values  y _j,y ^,y ^ / 2>y 5 / 2  at  heights 
z_|,Z|,z^y2»z5/2  ^as  f°r  above),  i.e. 


(■»*  =  ci'4-W'1 


(A.<r) 


where 


a+  ~ 


&+  —  ^/rx  ■+  y- 


(A.6£) 


(For  values  of  z^  other  than  70  km  (the  adopted  value),  (A. 66)  may  need  to 
be  replaced  by  an  alternative  formula).  Aq  is  a  relative  pressure 
increase  at  height  z^  that  may  be  associated  with  the  solar  cycle 
corresponding  to  (A. 36),  we  have 


a  =  aR„ 


(A‘7) 


where 


i  =  '  ~h'Cr  +  t  C^i) 


Gm*) 


^(zp.M  )f  the  diffusive  profile  number  density,  and  n  )» 

^  *  2 
the  N2  mixing  profile  number  density,  at  height  Z2  that  appear  in  (A. 59) 


_3 

are  MSIS-86  parameters  that  are  calculated  in  units  of  cm  by  the  MSIS-86 
formulation  for  the  required  time,  date,  latitude,  longitude,  solar  activity, 
etc.  In  units  of  mb,  we  have 

=  ,0<f  ^  ^  (A  .  6<?) 

where  T^)  is  the  MSIS-86  temperature  at  height  Z2  for  the  required  conditions 
and^^xs  the  reciprocal  of  Avogadro's  number  used  in  MSIS-86,  i.e. 

4J1  =  /.U*.  /0"27  (A.lty 

Hence  ^Ni=  6.02410  x  1026  (kmol)"1. 

2.  Density 

At  height  z  in  the  range  (ZpZ2>  density  p  is  calculated  from 

f  =  Z  M.  (V‘)  CA.71) 

1  1*6 

_3 

where  n(z,M^)  is  the  number  density  (m  )  at  height  z  of  gas  constituent 
of  molecular  weight  Mp  the  values  of  M^  being  taken  to  be  those  of 
MSIS-86,  i.e.  Mj  =  4,  M2  =  16,  M3  =  28,  M^  =  32,  M5  =  40,  M?  =  1  and 
Mg  =  14  corresponding  to  gases  He,  0,  ^  O2,  Ar,  H  and  N.  When 
different  values  of  Avogardro's  number  are  adopted  for  the  lower  and 
upper  models,  we  take  a  smooth  transition  given  by 

4«  =  i[  '4N,+'4t42  *  )tcini  cj]  (-  I  £  S  <  1 )  (A -72) 

where 

A  —  6.0225-7  >  lo 26  .</  -  (z.OZLt-lO  K  to  (4.73) 

N I  N 2 

in  the  present  calculation.  For  nCz.M^)  we  adopt  the  relations 

of  the  MSIS-86  formulation  and  have  (Appendix  A5) 


=  i°i[r'j(r1Mi)  +*y.lr>Mcy' !J  ‘Cji)C,A)  (A-74) 


A t  - 

wdtfy, -A)  =  nd  <■?*,”;)*-  frcr,)/ 1  by 


(A.75\) 

(A.7<0 
(A.  77) 


J"(J)  =  Luoi-uisj]^/^ 

2  x  t,s  -  .  - .  +  K-f7 

t>5  I  ' 

0=',^ 

.  /  -  ,  n  /  '  llT  .  _  _/  ) 

Y//(^'k)  C  n-v  1  rv-t-  3  n+5-  r>+7  ' 

ya/j)  =  /  -•/*+  |5+-  7  J-4 

where  o(^  =  -0.4  for  i  =  1  and  7  (He  and  H)  and  is  otherwise  zero. 

Cji  =  1  for  N2  (i  =  3)  and  =  1  for  He,  ^  and  Ar  (i  =  1,3,4  and  5). 

For  the  remaining  cases 

Cjity)- exp\rj:ll  i +■**[>  (j  =  ',2)  (’A-S3) 

l~bi  ~  R zi  C  ft.  9s) 

where 

R,;  -  (a.sO 

The  values  of  Si . ,  z..,  H..,  R0 .  are  those  of  MSIS-86  (Ref.  2, 

'■1’  ji  ji’  2i 

Table  2d)  and  for  i  =  2, 7  and  8  (0,  H  and  N) 

e,'-  =  A,i  ^.?7) 

whereas  for  i  =  1,4  and  5  (He,  09  and  Ar),  r0  is  chosen  so  that  these 


(A.  7  5?) 

(A.7<?) 
(A.  So) 

(A. Si) 

(A  .  §2) 


constituents  have  been  given  volume  fractions  of  air,  f^,  at  height  z^. 
The  required  values  are  given  by  (Appendix  A5) 


— 

[*,;  -  »: Afh  e ‘"It  W] 

(A.ss") 

N 

)  l(y3) 

(A.ai) 

!±c  =  p  pty )  / H<0 

(A.  7 '«) 

(A  .<5/) 

for  i  =  1,4  and  5,  where  we  take 

«*)-  t-.^0  -  >-000*13.  cmO 

In  order  to  maintain  continuity  in  the  mean  molecular  weight  at 
height  Zj ,  the  value  taken  for  is  such  that  the  mean  molecular  weight 
of  air  at  height  z^  for  the  gas  constituents  ^  C^,  Ar,  He  is  equal  to 
M^,  i.e.  such  that  28 fN  +  32fQ^  +40fAr  +  4fHe  =  Mz  <■  28.9644). 

The  value  obtained,  0.21023,  is  slightly  higher  than  the  actual  value, 
0.20948,  the  excess  being  largely  accounted  for  by  the  presence  of  CO2 
in  the  real  atmosphere  which  is  omitted  in  the  model. 

A  useful  check  on  computing  accuracy  is  obtained  by  evaluating 
for  i  =  3  (i.e.  N2)  and  comparing  its  value  with  the  adopted  value  of  zero. 


3.  Pressure 


Total  number  density  is  obtained  as 

yi  -  2?  wy  M,) 
Mean  molecular  weight  is  then 

M  f>  )  yu 

and  pressure  is 

j,v  =  c  £•/  ~r 


)  £m0 

(A}(k~et.Y')  {*.%) 
CA.?r; 


4.  Turbopause  height,  z^  of  the  i  th  gas  constituent 

In  the  MSIS-86  formulation,  turbopause  height  is  taken  to  be  the  height 

z^  at  which 

and  z^  is  a  fixed  parameter,  being  105  km  for  ^  C^,  Ar,  N  and  0,  100  km 
for  He  and  95  km  for  H.  In  the  present  formulation,  where  MSIS-86 
parameters  are  adopted  for  the  height  z^  which  (at  130  km)  exceeds  the 
heights  z^»  the  MSIS-86  values  for  z^ cannot  hold  unless  the  derived 
temperature  profile  in  (z^^)  *-s  t*ie  same  as  t*ie  ^^-86  profile, 
which  clearly  is  not  the  case.  However  as  deviations  between  the  two 
profiles  are  not  large,  the  values  of  z^  satisfying  (A. 96)  would  not  be 
expected  to  differ  greatly  from  the  MSIS-86  values.  In  the  computations 
that  have  been  undertaken  the  differences  have  been  less  than  1  km. 


Vie  calculate 


y*i  =  +  S*: 


(AH) 


where 


CM*) 


The  denominator  in  (A. 98)  is  calculated  from 


[vn)/s]  =jri.Jww-v'„y,(!)]r1* 

J  ^  is  obtained  from  (A. 98)  by  iteration  with  a  suitable  initial  value 
(corresponding  to  say  =  100  km). 


APPENDIX  A1 


FORMULATION  OF  W 


r  1 


We  define 


W~‘  =  A  +~5  (Al.l) 

where,  for  given  S,  N  and  a  (which  are  independent  of  f  and  5  ) 


sn 
S  Nf 


-B  -  £  ir5r" 

rs  t  r 


fA1.2) 

CA,.  3) 


^/n  are  chosen  so  that  for  all  values  of  jj 


r^i.o 


-i 


Hence 


v  =  ej<?. 


(Al.s) 


where 


e„-f  o-stj^j 

-I 

On  integrating  by  parts  it  may  be  shown  that 


CAf.6) 


C^+7)  ~  (^“0 


CA1.7) 


Hence  for  n  even 


/  =  C  €j  <2„-.  )...(£* /€„! 


_  1.3  (A  < . «  ) 

<^.ll  ,,  .  )o u  ) 


For  n  odd,  G^  =  0  by  direct  evaluation  and  hence  =  G^  =  ...  =0  and 

V„.  ~  0  L-n  arU')  C ftl?) 

^  are  chosen  so  that  for  given  . . . ,  H-j,  by  (Al.l)  to  (A1.4), 


1 
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(Al,  i<?) 


j'  8J5  =  f'w'jS  =  i, 

-I  ' 


b£0®  w"(-0^  j2  «<*#-•)-  A/7-0-  a 

*  J5  <*r  *'J?»  ~iJ'  S 

(Al.  n) 

B(0  =  w"'0  )  ~  ^3  j  <dJl^  ^  =  cLv)~'[  ij  -  Jl  oL*3(l  )  —  d\j  (i)  _ 

dS  O IS  ’  dS*  d  S  7 

(>4/.  /2) 

These  conditions  require  that 


where  b  = 


E^l  * 
$"*  - 


s  4  «  i 

-67]',  /-  ft  ...i7r 


and 


X  o  \  o  \  O  v7  1 

I  -I  1  -I  I  -I  / 

I  f  f  f  (  t  ( 

0  /  -2  3  'A  5“  -6 

O  /  2  3  4-  S"  6 

0  0  2  -6  /2  -5o  3o 

OO  2  (>  rz  2.  o  30 


C*M3) 


(AUA-) 


Fran  (A1.13),  we  obtain 


1  ~  ssi 


(At. 


where,  on  inverting  S 


-1 


5  = 


lO S  -9"  ->57  -«  «  -/  / 

^  -  */0  tjO  -42  -62  & 

l  Sts’  fo  11 

0  60  -  60  60  60/2  '12 
3/s  -3/r  -3/r  “wo  izo  -is  -/s' 

0  “/8  /8  -/?  “/g  -  6  6 

-/os-  /os'  /o*r  42  -42.  7  7 


~ 


(Al  ./6) 


APPENDIX  A2  CONDITIONS  AT  HEIGHT  zl 

By  definition 

w (A2-0 

ft  T  TCV) 

where  the  dependence  of  T  on  longitude  'X  is  now  indicated  and  Tref 
denotes  the  zonal  mean  of  T(  'X  ) . 

At  height  Zp  which  is  at  or  below  80  km,  we  have  R  =  Rp 
g  =  g^  (as  defined  by  (A. 6)  and  (A. 7))  and  from  Ref.  1,  Appendix  B 

HV- ^-2) 

A  T-  TOO-T^  -.i;crtOl->Xl)  +  T2«^C2>->T*)  w-V 
where  §  -  sin( latitude ) ,  Z  =  (Z-Z^/Z^,  Zg  =  48.75  km,  Z^  =  31.25  km 


31x1  cns’  Tl»  T2’ 

\r  x2 

are  tabulated  in  Ref.  1. 

Hence  from  (A2.1) 

to  (A2.3) 

W'1  =. 

where 

*r  = 

*»  = 

1  +  Ti 

v  ~ 

41r[Jr-(R,AT/Mli, 

(A  2.7) 

From  (A2.4) 

and  (A2.5) 

JW-1^ 

-8,  _  -fi,  J?t 

(A2.8) 

<4 

JD,  Jy 

and  from  (A2.6) 

and  (A2.7) 
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(A2.<)) 


Confining  (A2.8)  and  (A2.9)  with  the  use  of  (A2.6)  yields 


JUvJ‘l)/dy  -  (Aawo) 


where 

3>*  *  A'  U  ^2'"^ 


From  (A2.10)  and  (A2.ll) 


i^f'  =  -2j>a  4?,  +u  -iiA-  (*2-<2) 

a?,51  Jv  k  5  4 h,x 

by  <A2.5);  and  from  <A2.6)  and  (A2.9) 


dJ',/jy=  (A5.1S) 


vrtiile  from  (A2.5)  and  (A 2.7) 


dpIJy  -  +  (A2.  (*.) 

Hence  (A2.12)  to  (A2.14)  give,  on  using  (A2.ll) 


d'u-'/J.f  =  C3>,I>3+  2^') /(*,!>’')  CA2.IET') 

where 

3?3  -  iA-i-C-?, 


The  required  conditions  at  height  z^,  which  are  expressed  by  (Al.ll), 
become  by  (A2.4),  (A2.10),  (A2.15)  and  (A. 5) 
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-  -A,  /  X>,  f  A,  (A2.i 7) 

Ci>,y3>,x+Aa)^  cA2  '«) 

A  =  [  03,1),  +  2  J)2  )/ ,d,5)  *  A,1  £  <?*  • 


where  h  ,  qr  and  Dr  are  evaluated  at  z  =  and  the  terms 

Ar  -  [  rmy^\ 

are  introduced  to  enable  an  imposed  change  Aw  1  of  W  1  in  the  vicinity 
of  z  =  Zj  to  be  taken  into  account  due,  for  example,  to  a  dependence 
on  solar  activity.  Corresponding  to  an  increase  in  temperature  AT 

AW'1  =  -W'T”'AT  [fKTi.iO 


A  formulation  has  been  given  for  At  corresponding  to  a  change  in 
sunspot  number  ARn  (3),  by  which  (A2.20)  and  (A2.21)  yield 


Ar  ' 

Kef) 


Kff)  ir  /\R„ 


n  , 

c<rS  4> 


g  _  £ 


r  m 

ir-'  [|t,W’2 


-i 


^2. 


where 


CA2-2*z) 

(A2.2?) 


2.  =  (A2-2S) 
The  quantities  <*,  a,  p,q  and  n  are  calculated  from 

s  =  e,  i  0*.  °*"n)  (A2/2i) 

|  =  taJhik-fln)  ca[ltCtj-Olls2  Sl  [*7-2 1) 

where  t^  is  the  day  number  in  the  year  and  =  0.60,  =  66.63  km, 

=  12.14  km.  Pi  =  7.19  x  10  ^  km  ^  q^  =  4.80  x  10  ^  km  1 ,  =  6.0, 

<^2  =  0,  a2  =  -3.09  km,  ^  =  0,  P2  =  -0.51  x  10  ^  km  q2  =  -2.16  x  10  ^ 
km"1,  n2  =  -5.0. 
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APPENDIX  A3  CONDITIONS  AT  HEIGHT  z: 


At  height  Z2  (^ 90  km),  MSIS-86  is  valid  and  we  have  R  =  R^,  g  =  ^  (as 
defined  by  (A. 6)  and  (A. 10))  and  (A.l)  becomes 

(A5-° 


By  logarithmic  differentiation 


c w-r 

=  -2^ 

-  T"  J-T/Jy 

£A».») 

where 

<L  - 

■  '/  (R 

(■+%) 

(A3.  3) 

A  further  differentiation  gives,  by 

(A3. 3), 

i’W1 

-  (w- 

y\%T ' 

2  4J)z__i 
l  J  | 

(A3.*) 

Hence  from  (A3 

.2)  and 

(A3. 4) 

(W 

dL’WT’ 

(0 

d  +  4  ^  ^ 

1 

_i  d!Tf  _  j.  o(ar 

T  /  T  dj* 

(A3.?) 

2  2 

Different  expressions  hold  for  dT/dz,  d  T/dz  according  to  whether  z  is 
greater  or  less  than  z  =  117.2  km. 

For  z  ,  we  have  (2) 

Z  d 

712,)=  (A,'0 

Hr, ft)  ~  ^'VA)t^|,7^)/C^  +  5r)  (7*V7) 

where  T^  is  exospheric  temperature,  is  temperature  at  height  z^ 

(=  120  km)  and  R  =  6356.77  km.  c  is  a  constant  that  is  related  to 
P 
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the  temperature  gradient,  T^' ,  at  z  =  z^  by  =  'T^/(Tee-T A  )• 


From  (A3. 6),  we  have 

4g"V-?-> 

(A3.*) 

(AS?) 

where 

=  a*  <*!  Uy 

(A3.  io) 

and  from  (A3. 7),  by  (A. 10)  and  (A3. 3), 

j?Uy-=  (Rr+?u)7(V})'= 

(A3. if) 

yoi  - 

(A3.)2) 

Hence  if 

*  =-«-3Ci--^/-r) 

(A3.I3) 

Yj  —  -U.  (2v  +  5"^ 

(A3.  /a) 

(A3. 2)  and  (A3. 5)  give  by  (A3. 8)  to  (A3. 13) 

(vT*)'1  cl  W'7<^  =  -(^c{4--u) 

(A?./?-) 

=  Gcf+Gv  X-r-T) 

(A  3.  >0 

where  by  (A3. 10)  and  (A3. 11) 

«^=*  <r  3,<»)/?^yAN) 

(A3. 17) 

For  z  ,  we  have  (3) 

Z  « 

r"'  =  ro'V  rc  (A3./8) 

where 
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*  -  -  r t<-r,u)- siv.,\y 


so*,*.)  =  c>-vKVk)/0V>) 


(A3.2o) 


zq  is  mesopause  height  as  determined  by  the  MSIS-86  formulation.  By 
(A3. 18)  we  have 

Jl4T  — +5rj>*V£  &*•*'> 

Tdy  d*  J" 

2  ( 1  <Ll  ~\  -  J- 4Ll  -  T  i-1  -  53cT(Tg<-2TaT+3't)3+)—  -  [jh(—  jl 

2WdjJ  Ta^  av"  ~  c  '<W  M 

+2T(^6T0*+)T^x)^\’'  (A*.2») 


(A3. 2)  and  (A3. 5)  can  then  be  written,  by  (A3. 19)  to  (A3. 22)  and  (A3. 3), 


as  (A3. 15)  and  (A3. 16),  where 


rj  =  JT(T6+;rcxt 


(A  3.  .23) 
(A3. 74; 


and  by  (A3. 19),  (A3. 20)  and  (A. 10) 

0(31  __  _  j  ) 

~  Jaffc) 


(A3.25) 


The  required  conditions  at  height  z2,  which  are  expressed  by  (A1.12), 
becotre  by  (A3.1),  (A3. 15)  and  (A3. 16),  on  changing  units  from  m  to  kra  in 
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A  =  loJ  Mh  /^T(v.)  (X.") 

(AwT'J 

Jlj  —  (  6  +  £  W.  o(  ■+•  ^)  ^ 


6*3.26; 

(A-3.  27; 
(>*\20 


where  d,  u  and 


are  evaluated  at  z  =  ^2' 


APPENDIX  A4 


THE  N2  PRESSURE  RATIO  CONDITION 


We  adopt  the  MSIS-86  representation  for  the  number  density  profile  of 
a  gas  constituent,  n(z,M),  as  a  meld  of  a  diffusive  profile  n^(z,M) 
and  a  mixing  profile  nm(z,M).  For  molecular  nitrogen,  M  =  and 
the  representation  at  z  =  z^  is 


~  inatooM.<I) 


(A  4-0 


where 


A*  = 


(A4-23 


and  =  28,  Mq  =  28.95  kg(kmol)  On  multiplying  (A4.1)  by 

jjR-^/4^)  T(z1)jAN2>  R1»^N1  1)6:5118  defined  by  (A. 6)  and  (A.  73),  we 
obtain  a  relation  in  terms  of  the  N2  pressure  p(Zpf^  )  at  height 

A)*"*  -  Mir, ,  0*4-5) 


where 


kCy,M.)  =  *,(*.  «;)T(») 


(M-.iL) 


with  ^ 


By  integration  of  the  hydrostatic  equation,  the  diffusive  and 
mixing  pressures  are 

-M;I  -|0(;  /  . 

*MM)  -  (Afc.7j 


where 


ito.to)  =  Jj- 


On  dividing  (A4.3)  by  |pm^z2»^  )]^2  we  obtain,  using  (A4.4)  to 


(A4.7)  with  =  0 


(M0/M*)X,  A,  _X  \A„l 

'  ~  (VN.  ^  +  V 


(A  A  ■)) 


where 


X  =  i  (y^) 

At,  = 

VN_  ~  A*  fra,  M^v)/  ^^5,  ) 


(A4.ic) 

A  4-  -  i1 ) 
(AA  .  iz) 


fN  is  the  ^  volume  fraction  of  air  at  height  z ^  and  p(z^)  is  air  pressure 
at  height  z^,  which  is  evaluated  as 

[*<*,)=  (a^iO 

where  pre^  is  the  zonal  mean  value  for  a  reference  value  of  solar  activity 
(sunspot  number)  and  Dq,  Aq  are  relative  incremental  increases  associated 
with  longitude  and  sunspot  number  respectively.  We  have 


3> 


-  [|>M-/>.,f)/f.r^  ’  (A4'4;> 


where  p^,  P2»  ^pi  ’  are  tatiu’ a' t-vi  :n  Rfc,f-  1-  Corresponding  to 


(A2.22) 


=  Kit)^  i,  A? 


(^A  U-.iS) 


where 


i„  -  [  ^*0^-4-,)] 


(A  4  .  16,) 


The  required  pressure  ratio  condition,  which  is  expressed  by  (Al.10), 


is 


A  =-  (l\/Mj  X  /», 


(AA.'7) 


by  (A.l) ,  (A. 5),  (AA.8)  and  (A4.10). 
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APPENDIX  A5 


NUMBER  DENSITIES 


For  a  gas  constituent  of  molecular  weight  M^,  the  MSIS-86  representation 
of  number  density  n(z,MJ  introduces  factors  C^(z),  (^(z)  to 
account  for  effects  of  chemistry  and  flow: 


Cjc l  ?r)  ~  \  rji  /  [>+  £>  -y-ji  1}  (j=^2>  CA5: ') 

=  1  for  N2  only  (i  =  3)  and  C2i  =  1  for  all  gases  except  0,  H  and  N 

(i  =  2,  7  and  8).  Hie  values  adopted  here  for  z..,  H..  (i  =  1  2)  and 

J1  J1 

are  those  presented  in  MSIS-86  for  all  gases.  So  also  are  the  values 
adopted  here  for  r^  (i  =  2,  7  and  8)  as  these  gases  (0,  H  and  N)  are  not 
required  to  match  significant  volume  fractions  of  air  at  height  z^ 

(=  70  km).  The  other  gases  (He,  O2  and  Ar)  are  however  required  to  match 
significant  volume  fractions  f^  (i  =  1,4  and  5)  at  height  z^  and  this  is 
achieved  by  means  of  the  relation 


r,i=  (Am) 


valid  for  where  c  and  £  are  defined  by  (A. 4)  and  (A. 3). 

R-^  depends  on  f ^  and  is  derived  as  shown  below.  Values  of  R-^  computed 
in  this  way  are  found  to  lie  within  about  one  per  cent  of  R^.  The 
adjustment  using  (A5.2)  to  provide  a  smooth  transition  between  lower  and 
upper  models  (for  He,  O2  and  Ar  concentrations)  is  therefore  of  negligible 
physical  consequence. 

For  the  calculation  of  n(z,M^)  we  have  (2) 


At' 

n<r,Mi) =  Lw,  ir,  ;]  ‘  c„. ij)  c„<»i 
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A;  - 


(p~-3)  C«r.r) 


McJ(S)  %  _  . 
TM/T(v) 


Lc~-")  (AS-.4) 


where  o<^  is  the  thermal  diffusion  coefficient  and  by  (A4.8),(A.l)  and 


(A. 5) 


oth  =  i(r,  *,)  =  i«o  -</«)>.,  I  rtti 


where 


0(5)  =  j^w"  JS 


(AS--7) 


(A^-O 


By  (Al.l)  to  (A1.3) 


UC5)  -  X  X  2  ^ 

-a*i  n  =< 

+jcs)-STa-rnr-<)JS  ^w) 


^[^T- 


3J_Z 

v»+3 


WJ)  =  «  -  ^  * 


1  5  * 

5"  *> 


n+4 


—  r6 
7  •> 


Q'VS'-m^ 

(frr.ia) 


To  evaluate  (i  =  1,4  or  5)  we  note  that  C2^  =1  (i  =  1,4  or  5) 
and  put  z  =  z-^  in  (A5.3)  to  obtain,  on  multiplying  through  by 


4* 

V 

y 


I 

| 

ii 

i 


E 
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(Rj/>rNj)  T(Zj)/pm(z2,M^)  and  using  (A4.4)  to  (A4.7)  and  (A5.7), 


^  ;  lv-  * 


f^*TTU«  <* 


where 


Hi  =  fit*)/  Mfr) 

v.  =• 

By  (A5.1)  and  (A5.13),  we  may  solve  for  r^  (=  R^),  where 

R,-  =  £  I  i-  ^[CY'-lT'cl/h)' /]  ^  JL  c(£  (y.) 


J.  <ww  -  A.  „  -  4?Ajh,w-°( 5-*f  ]';  t-^T  1 


APPENDIX  B 


METHOD  OF  DETERMINATION  OF  a 


sn 


We  write  (A1.2)  as 


S  *J 

2  Z  ■=  A 


--4  “»  ■«  - ' 


C»-0 


where 


CKO  - 


-4-1 


and  determine  agn  by  the  method  of  weighted  least- squares  from  sets  of 
values  C^s.n)  (s  =  n  =  1,...,N),  for  k  =  1,...,K.  The 

weighting  is  based  on  estimated  standard  deviations  of  A^,  irA^,  which 
are  obtained  as  described  below. 

The  K  sets  of  values  are  taken  at  grid  points  located  on  a  meridional 
cross-section  at  each  5  km  height  interval  from  to  Z2  (excluding  the 
end  points  where  11  and  no  contribution  is  made  to  the  determination  of 
a  )  and  at  each  10°  latitude  from  pole  to  pole.  Then 

Oil 


K  —  '1  [(Y.-V,)/S-  -  ‘J  (.8.5) 

For  the  adopted  values  z^  =  70  km,  Z2  =  130  km,  we  have  K  =  209. 

For  each  month,  are  determined  as  the  weighted  average  of 
differences 


A  -  Id T  -  B 


where  T  is  observed  and  B  is  calculated  for  the  same  geophysical 
conditions  as  the  observation.  When  calculating  B,  the  terms  which 
introduce  an  observation's  longitude  and  incremental  sunspot  number. 


&Rn. 
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into  the  lower  model  (as  defined  in  Appendices  A2  and  A4)  were  found 

to  be  small  and  ineffective  (in  the  a  determinations)  and  were  therefore 

sn 

omitted  in  obtaining  the  agn  values  presented  below.  Longitude  and 
solar  activity  dependences  were  however  retained  in  the  upper  model 
(Appendices  A3  and  A4)  in  the  calculation  of  B.  Solar  activity  is 
expressed  by  the  F^q  j  ,  F^q  j  and  A^  indices  which  were  assigned  to  the 
temperature  data  as  follows: 

(i)  Time-averaged  data 

Incoherent  scatter  monthly  means:  Millstone  Hill  Ap  =  7,  F^g^=120; 

St.  Santin  Ap  =  7,  Fio  7=70>  170  (i*e.  three  sets  of  monthly 

means);  Arecibo  (1970-  75)  Ap  =  10,  F^g  ^=108  units. 

Rocket  monthly  means:  Kwajalein  (1976-  78)  Ap  =  10,  F^g  y=95;  other 
sites  (none  of  which  provide  data  above  100  km)  A^  =  10,  F^g  -j  =  120  units. 
Fjg  -j  is  the  3-month  average  of  F10.  -J  and  was  taken  equal  to  F10.7 
for  time-averaged  data.  A Rn  (when  not  set  equal  to  zero)  is  calculated 

£ran  F10-7  by 

fiR^  =  I. os  FK1  -  61  -  R„0  CB-S") 

where  RnQ  is  the  reference  sunspot  number  appropriate  to  the  lower  model 
and  taken  as  65. 

(ii)  Single  rocket  profiles 

From  the  date  of  launch,  R^  is  calculated  by  linear  interpolation  of  a 
table  of  monthly  mean  sunspot  numbers,  then 

F,o.7  =  (R l-os 

Likewise  F^q  -j  is  obtained  from  the  3-monthly  sunspot  number.  Ap  may  be 
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specified  for  any  launch  or  otherwise  set  equal  to  10.  Few  profiles 
extend  above  100  km  at  high  latitudes  (see  Table  3)  where  dependence  on 


Ap  may  be  important.  After  examining  A^  values  for  a  sample  of  such 


launches,  the  approximation  =  10  for  all  launches  appeared  to  be  justified. 


Dependence  on  local  solar  time  in  the  upper  model  is  included,  when 
known,  as  in  the  case  of  single  temperature  profiles.  For  monthly  mean 
temperature  data,  B  is  calculated  for  the  diurnal  mean  upper  model  with 
the  exception  of  incoherent  scatter  data  which  are  obtained  during  daytime 
hours  only.  For  I.S.  data,  B  is  evaluated  for  each  hour  between  0800 
and  1600  hours  and  the  average  is  taken  of  these  11  values. 

The  weighted  average  of  the  values  of  A  calculated  frcxn  (B.4)  for 
the  k  th  grid  point  is  based  on  those  temperatures  at  the  height  of  the 
k  th  grid  point  that  lie  within  15°  latitude  of  the  k  th  grid  point  and 
whose  date  of  observation  lies  within  1.5  months  of  the  middle  of  the 
given  month.  The  weighting  factor  used  in  the  averaging  process  is 


e-x 


p  \  -  i  L£W=/'°)  + 


t*.7) 


where  is  the  latitude  displacement  from  the  grid  point  (in  degrees)  and 
Am  is  the  date  displacement  (in  months),  A-^/10,  Am  being  less  than  1.5. 
The  most  extreme  data  point  to  be  included  in  the  average  is  therefore  at 


=  15°,  Am  =  1.5  and  has  a  weight  of  e  (=  0.11)  compared  with  a 


data  point  at  the  grid  point  (  tSj>  =  m  =  0)  for  which  the  weight  is  unity. 

Temperature  data  are  available  either  as  single  observations  or  as 
monthly  mean  values,  but  the  relative  accuracies  of  the  various  types  of 
data  are  invariably  unknown.  The  procedure  adopted  for  assigning 
relative  weights  has  been  to  weight  all  single  profiles  equally  and  to 
give  all  monthly  means  3  times  the  weight  of  a  single  profile.  A  factor 
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of  10  (instead  of  3)  was  initially  adopted  but  this  resulted  in  data 
for  single  profiles  being  almost  completely  discounted  by  that  of  the 
monthly  means  in  cases  where  both  were  combined.  The  value  of  3  was 
therefore  adopted  and  seems  to  allow  a  more  realistic  utilisation  of 
single  profile  data  in  the  analysis.  The  final  determinations  of  a£ 
are  not  sensitive  to  the  weights  adopted  as  data  for  single  profiles  and 
monthly  means  are  from  different  locations  and  tend  not  to  combine. 

The  introduction  of  the  weighting  factor  (B. 7)  results  in  fractional 


sn 


values  for  the  weighted  number  of  observations,  N  ,  and  the  minimum 

w 


value  for  Nw  below  which  the  available  temperature  data  are  considered 
insufficient  for  the  determination  of  A^  is  taken  to  be  1.1.  In  that 
case  is  taken  to  be  the  value  of  A^  for  the  same  latitude  in  the 
opposite  hemisphere  with  the  month  shifted  by  six  months,  if  such  a 
determination  is  possible  with  the  data  available.  By  this  procedure, 
it  was  found  that  A^  could  be  determined  at  all  grid  points  except  for 
a  number  at  high  latitude  (60°  or  more)  and  above  100  km.  For  these 
grid  points  A^  is  determined  from  (B.4)  with  T  equal  to  its  MSIS-86 
value. 

At  the  k  th  grid  point,  crA^,  the  standard  deviation  of  the 
weighted  average  A^,  is  estimated  from  the  standard  deviation,  0“,  of 
the  distribution  of  values  of  A  in  the  average.  The  relation 
adopted  is 

(*•*) 


=  <r/  (tju-  uy 


in  which  cr  is  also  an  estimated  value.  In  those  cases  where  A^  is 
found  from  MSIS-86  temperatures,  <yA^  is  calculated  as 
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„-A*  <r(Ma/KT)  =(Mti1/RTJ(ffT/Tj  (e.<j) 

where  Tc  is  the  temperature  value  from  CIRA  1972,  Part  1  and 

je^(cT)  —  o~  (y -47. S’)/ 4-2. S'  -+  £ 

where  a  =  0.519,  b  =  0.617,  (being  values  based  on  CIRA  1972  temperature 
distributions  (4)).  (B.10)  gives  <tT  =  4.4,  10.3,  24.0  K  at  70,  100, 

130  tan  respectively. 

In  those  cases  where  <rA^  can  be  found  from  (B.7),  it  was  desirable 

to  avoid  unreasonably  small  values  of  <TA^  which  may  at  times  arise 

with  selective  distributions  of  data.  Hence  is  compared  with  sA 

and  if  found  to  be  less  than  <rA/3,  we  take  flrA^  =  <r A. 

The  procedure  for  determining  S  and  N  in  (B.l)  is  to  take  them  as 

small  as  possible  and  consistent  with  a  satisfactory  least-squares  fit  to 

the  values  of  A^.  The  values  chosen  are  S  =  7,  N  =  2. 

Sets  of  values  of  a  determined  for  each  month  for  Case  2  of  Section 

sn 

4.2  are  listed  in  Table  B.l.  Case  2  omits  all  data  at  and  above  105  tan  in 

order  to  get  an  improved  fit  to  data  at  75-  100  tan.  The  number  of  grid 

points  of  the  height-latitude  cross-section  used  for  the  fit  is  144,  those 
omitted  being  at  the  five  heights  105,  110,..,  125  km  and  at  the  13  latitudes 
60S,  50S,..,  50N,  60  N.  Grid  points  above  100  km  at  latitudes  70,  80, 

90°  N  and  S  are  retained  so  that  the  polynomial  fit  gives  reasonable  values 
at  high  latitudes  above  100  km. 

Case  3  which  omits  all  data  below  100  tan  in  order  to  get  an  improved  fit 
at  100-  125  km  utilizes  grid  points  at  70,  80  and  90° N  and  S  at  all  heights 
(75-  125  tan)  so  that  the  polynomial  fit  gives  reasonable  values  at  high 
latitudes.  The  number  of  grid  points  used  for  the  fit  is  133,  134  or  139 

according  to  the  month. 


Table  B.l 


Coefficients  a  calculated  for  Case  2. 
sn 


COEFFICIENTS  all 


JA'JUAftY 
all# ...#a7l  = 
*12#..., *7 2  = 

FEllfJUARY 
a 1  1 # • • • # *7 1  s 
a12#«.«# a? 2  s 

KAHCtl 

«»1 1  # . . . #  a7 1  = 
*12#..,#o72  = 

APRIL 

*1  1 . *7  |  = 

*12#... #*72  = 

MAY 

all#. ..##7  I  = 
a 1 2  # . . . #  *72  ” 

Jl'HE 

all#...#  <7 1  — 
*12#. ..,*72  = 

JULY 

all#. ••#*71  = 

*12#. ..,*72  = 

AUGUST 

a 1  I # . . . # *7  1  - 
ol  2 . . i72  = 

scpterucm 

all#.  ..,*71  = 
al 2# ... #*72  s 

or  rcuEM 

all#  ...#*7I  s 
*12#  ...#*72  s 

uC^CmBEP 
all#.  ..#*71  s 
a  12#. ..#*72  = 

oectME* 
all#..., a?  I  = 
*12#. ..,*72  s 


0.4502370*01 

•0.52?20o£4C1 


0.4716V*EtC1 

•0.3of5270401 


n.43?on9F+ni 

•G.4796iaF»nu 


0.8944430.401 

o.i2<7«>4*;E»ni 


n.6o476SFtm 

O.Ao°U«oEtCU 


0.61723„E*C1 

o.iojaoiotoi 


0.54?099E*01 
•0.1  c»7ol  9E  *01 


n.Sj704i)Ft01 

•0.1u«923EtG1 


0.4990^oCfG1 
0 .8U4<«7j«E  tOO 


0.54C66uEt01 
0.?<j9o2VE  tOl 


0.41 3  78VE  *01 

n.i4ftr>22Etni 


0.3l530iCt01 

0.SU27dE*01 


-0.4J1295Et01 

0  -  4  JS  1  s..c  *0  I 


-o.l  I2.J49E t02 
0.1'J8*0jf  tOc 


-0 .4  J8  J3  jF  +0  | 
0.4*2 J6JF  tO  < 


0.1 90*7 JF tO J 
C.  124 >4  VC t02 


-0.4'*7  jCJEtOj 

0.1 #778jEt02 


C .  1*8*,5lF  t0  | 
C.1  lAjOXtOj 


C.3j499JEtO| 

-0.526j9uEt0l 


0 . 2J0 J8  2f  tOl 
-0.7  j956,E»0i 


-0.1 J7jOaEtO/ 

0.1  l7j7.,Etm 


-C. 47690  .EtOl 
-0.4J679  |E  »0l 


-C . 4  J477*#C #0l 
-0.1 20*73 7F  tOj 


-0 . 1  «i?I»2 1C  tO| 
-0.721  J1  lO.O.j 


0.410995E401 

0.79«57bEtO1 


0.1iS038Et02 

0.3*6494E*ni 


0.127738E402 

0.8.'5«.62C401 


C.828l47r.#ni 

0.3  i2U6lE«02 


■fl  .8  70562C  400 

•0.2  I410JC  #0u 


C.1 17084E#01 
0.591o46£f0l 


•0.2.#0554C401 
0.1 *6J29E*02 


0.2  I3179E401 
0.1a7024E402 


.#•71#*  *12,. 


0.24837lEt02 

0.2*22390*02 


0 .679V79F 1  CO 
0.1,,0JcCFiC2 


0.7140530  4OI 
0.4 102660  402 


0.1I#658lEt02 

0.229324F #02 


•0.226V9UE  4  02 

0. 2u7j99E  4  02 


0.1 J3521F#C2 
0. 1 J8*64F  t02 


#a7  2 


0.8714/40  tOU 
0.66249lEt02 


-n.131 20ir  #02 

-0.1b870&r.  #02 


0 ,57229©E  4OI 
0.743114E402 


-0.2?43l2EtC1 
0. 738/5tlE  »02 


0.148o96E402 

0.105370C#0J 


0.4046530*01 
0 . 1 27&8VF  *03 


-0.455yi4Ft02 
0.135V3VF  *03 


-0 . 564U56E  *02 
0 . 2  7533oF  »02 


0.4187SSF #02 
0.*>6738u0*02 


-0.528824.0  #02 
0.3i42S<.E#02 


-O.l74i#4JEt02 

0.227.,A9Et02 


-0.3u?vnjC402 

0.2o327VE»02 


-0.1649° IE #02 
-0.1 JEt02 


-0.1 J?*22r*0l 

n.3  i?;72Ct0i 


0.4 33j7JEtO  I 
0.2l8/12Et02 


0.92233JE401 

n.5iJ6j7/Et01 


0.1 u?287E  t02 
-0.2u8^9/E»02 


0.1##1  J2VE*02 
-0.H8l62E#02 


0.9u9 j29E#0  I 
0.2285510402 


-n.5J0^siE#n» 

0 .  8!#?*2*0  tO  I 


-0.11 6-#°  JO  #02 
0.714J87E»0I 


-0 . 753  70 JO *01 

0.21°*890*02 


0.2J5646E402 
0.1 7218104 02 


0.2*25°70402 
0.9  75U34E#01 


0.279561E402 

0.971447E401 


0.2*27810402 
0.2  17S?5E#02 


-0 . 8 } 6 v8 3F ♦ f 1 
-0.2  I3‘.5UE#02 


fl.11«27«f #0? 
-0.22°204F  tCI 


0 . 2 J0482E  #02 
-0.5247530  tOI 


0 . 1  *#0  J830  1 02 
-0 . 1 o7 1 33F  *02 


0.427345E#C2  -0.2545*20402 
•0.4  25J96E  402  0.86067lETC2 


0.215197F«C2  -0.1 58VJ>5E  402 
0.417132F  4OI  0.701UO!>Et02 


-  59  - 


0.1 JS134E402 
0 ,7?564<.E  402 


O.?«64l3Et01 

r.9/S830Et02 


0.SJ232UE401 
0.9 iSo7oE402 


0.3J4U4AF401 

0.2-,°?1oE«02 


0.1  J24P.4F  #02 

C.8J4»ME#02 


n.°2°2‘3UE401 

C.7/°U81C402 


0.218331E402 

0.®o8af3F#02 


0.1..4VP1F#C2 
0.1 175800  *03 


0.1 372200*02 
0.1*12120*03 


0.2222470402 

0.42182OC402 


0.2  •#7/830402 

0.9*59730*02 


0.1901810.02 

0.59865 70 *02 


APPENDIX  C  AVERAGE  TEMPERATURE  DEVIATIONS  FROM  MODEL  VALUES 

AND  THEIR  MEAN  WITH  RESPECT  TO  DIFFERENT  SITES 

Comparisons  between  observed  temperatures,  T,  and  model  values  Tm  are 
analysed  in  terms  of  the  differences 

=  r  '  (C./) 

at  each  5  km  height  interval.  For  the  i  th  site  we  obtain  the  average 
deviation  at  any  height  as 

where  N^  is  the  number  of  observations  available  within  a  selected  group  of 
months.  Four  such  groups  are  considered,  namely  (i)  the  three  'winter' 
months  (DJF  in  the  N  hemisphere);  (ii)  the  six  'equinox'  months  (MAMSON); 
(iii)  the  three  'sunnier'  months  (JJA  in  the  N  hemisphere);  and  (iv)  all 
months.  The  standard  deviation  of  x.  was  estimated  in  each  case  as 

<r*t  =  <^) 

for  2. 

The  mean  of  with  respect  to  those  sites  for  which  a  value  of  x^ 
could  be  determined  was  obtained  as 

;=  z  w; y  «?■<*) 

•r 

where 

-urL  =  d*  s-*:)"*  (p’5') 

and  N  is  the  nunber  of  such  sites, 
s 

The  standard  deviation  of  x  was  estimated  as 

•■x  =  «-p  /CH,-*)4  0=-*) 

where  is  the  standard  deviation  of  the  distribution  of  the  values,  x^, 
about  being  estimated  (for  N  ^  2)  from 

b 


APPENDIX  D 


CODING  OF  THE  MODEL  FORMULATION 


In  the  course  of  this  project  Fortran  programs  have  been  written  and 
tested  for  each  stage  of  the  model  formulation  and  the  work  has  proceeded 
to  completion  with  the  confidence  that  the  formulation  was  computationally 
practical  and  efficient. 

Over  50  programs  have  been  developed,  some  of  which  are  subprograms 
to  the  main  calculation  while  others  have  served  only  a  transient 
purpose.  The  latter  category  includes  test  programs  (i.e.  main  programs 
for  temporary  use  in  the  development  and  testing  of  the  subprograms)  and 
programs  for  one-off  calculations  such  as  the  inversion  of  matrix  (A1.14) 
and  the  calculation  of  Millstone  Hill  and  St.  Santin  I.S.  temperatures 
from  formulae  in  the  references  quoted  in  this  report. 

The  main  programs  will  be  described  briefly  in  this  Appendix 
to  indicate  sane  of  the  chief  computational  stages  of  the  work. 
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Details  of  the  subprograms  and  input  files  that  are  utilised  ini 
conjunction  with  these  main  programs  are  given  in  Appendix  E.  The 
parties .ar  way  in  which  the  computation  breaks  down  into  different 
subprograms  reflects  the  stage-by-stage  development  of  the  method.  No 
retrospective  consideration  has  been  given  at  the  present  time  into 
restructuring  the  coding.  Computational  efficiency  nevertheless 
appears  to  be  quite  good.  All  computations  have  been  carried  out 
interactively  on  EUCLID,  (the  GEC  machine  at  University  College  London). 


1 .  Determination  of  a_n  (Appendix  B) 

TDIFAVP4  (Temperature  Differences  Averaged  Program  4)  calculates  A  from 
(B.4)  for  each  observed  temperature  at  5  km  intervals  and  takes 
weighted  averages  to  obtain  A^  at  each  grid  point  of  the  height- 
latitudinal  cross-section.  cA^  is  calculated  by  (B.8).  When 
temperature  data  are  lacking  A^  is  calculated  from  (B.4)  for 
heights  above  90  km  using  the  temperature  value  of  MSIS-86  and 
tf-A^  is  used  as  a  flag  and  set  equal  to  zero. 

(Note:  the  subprograms  listed  below  need  to  be  added  to  the  main 
program  in  all  cases.  Subprograms  that  are  already  contained 
in  a  main  program  are  not  listed,  being  usually  just  short 
routines) . 

Subprograms:  P0LYZ12S,  BCZ12S,  BCATZ1S,  BCATZ2S ,  EARTHS,  GTS5S, 
MSISINS,  SCATZ1S. 

Input  files:  TPW12Q,  COEFFDQ,  RZDATA,  THADATA,  POLYTDD,  DUMMYAXD. 
Output  files:  TDIFWj),  TDIFtN^Q,  TD(N1)SD. 
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WADJP6 


WADJP7 


(N^)  stands  for  an  integer  (the  run  number  and  is  used  to 
identify  the  outputs  corresponding  to  different  inputs). 

TDIF(N^)  contains  A^  and  TDIF(N^)Q  is  an  inspection  file 
for  output  from  the  various  stages  of  the  computation;  and 
TD(Nj)SD  contains  the  average  temperature  deviations  from  the 
fitted  profile,  x^,  and  their  standard  deviations, <rx^,  for 
the  i  th  site  at  each  5  km  height  interval,  the  average  being  taken 
four  times  with  respect  to  data  falling  in  each  of  the  groups 
of  months  (i)  DJF,  (ii)  MAMSON,  (iii)  JJA  and  (iv)  all  months. 

(W  1  Adjustment  Program  6)  obtains  a  by  weighted  least-squares 
solution  of  (B.l)  for  each  month  of  the  year.  Before  doing  this, 
it  processes  A^  and  5~A^  for  each  grid  point  of  the  height- 
latitudinal  cross-section  as  described  in  Appendix  B:  (1)  if 
A^,  <rA^  have  been  calculated  from  observed  temperature  data  they 
remain  unchanged;  (2)  if  <rA^  =  0  (i.e.  data  are  lacking),  A^, 

<rA^  are  taken  as  the  values  for  the  grid  point  at  the  same 
lati  nde  in  the  opposite  hemisphere  with  a  6  month  shift  of  date, 
unles.  this  also  has  trA^  =  0  (i.e.  data  are  lacking  there  as 
well)  when,  at  the  original  grid  point,  A^  is  left  unchanged 
(having  been  based  on  an  MSIS-86  temperature)  and  ^A^.  is  set 
equal  to  <TA,  the  value  calculated  from  (B.9);  and  (3)  any 
crA^  which  is  less  than  <tA/3  is  considered  improbably  small  and 
is  replaced  by  <r A.  WADJP6  is  for  Case  1  of  Section  4.1. 

Subprogram:  NAG  library  routine  F04ARF. 

Input  file:  TDIF(N^) 

Output  files:  OUIKN(N1),  WAIXHN^D. 

OUTRN(N^)  is  an  inspection  file  and  WADJ(N^)D  contains  agn. 
is  WADJP6  modified  to  deal  with  Cases  2  and  3  of  Sections  4.2  and  4.3. 
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2. 


t»«HH 


iarison  of  observations  with  derived  model 


TDIFAVP4  re-run.  TDIFAVP4  is  run  as  above  with  the  following  changes: 
Input  files:  as  before  with  DUMMYAXD  replaced  by  WADJ(N^ )D. 
Output  files:  TDIF(N2),  TDIF(N2)Q,  TD(N2)SD. 

Observed  temperature  data  are  now  compared  with  the  final  model, 

calculated  from  (A.l),  (A. 12)  and  (A. 13)  using  the  determinations 

of  a  which  are  read  in  from  WADJ(N-,)D,  whereas  in  the  first 
sn  1 

use  of  TDIFAVP4  above  a  were  unknown  and  were  read  in  as 

sn 

zeroes  from  DUMMYAXD.  The  purpose  of  this  computation  is  to 
obtain  TD(N2)SD  for  x^,  <rk^  for  use  with  AVTDSDP8. 

AVTDSDP8  (Average  of  mean  Temperature  Deviations  and  its  Standard 

Deviation  Program  £5).  For  each  site,  the  mean  temperature  dev¬ 
iations  from  the  model  with  respect  to  data  taken  in  each  of 
four  groups  of  months,  namely  DJF,  MAMSON,  JJA  and  all  year,  and 
estimates  of  their  standard  deviations  trx^  are  read  from  TD(N2)SD. 
The  sites  are  put  into  three  latitude  groups,  0-30,  30-50,  and 
50-90°N  or  S  and  their  summer,  equinox,  winter  and  all-year 
mean  deviations  are  averaged  and  the  standard  deviations  of  these 
averages  are  obtained  from  equations  (C.4)  to  (C.7).  Results 
are  shown  in  Tables  4-10  for  the  all-year  mean  deviations. 
Subprograms :  none 
Input  file:  TD(N2)SD. 

Output  file:  AV(N2)Q. 

3.  Tabulation  of  atmospheric  properties  at  heights  above  18  km 
TVALP  (Temperature  Value  Program)  was  written  to  develop  a  method  for 
combining  the  outputs  of  the  three  models  for  the  regions  18-70, 


70-130  and  above  130  km  for  given  ranges  of  geophysical 
parameters  (height,  date,  location,  solar  activity  etc.).  The 
program  was  first  written  for  temperature  only  and  later 
extended  to  provide  composition,  density  and  pressure.  Five 
subprograms  are  included  in  this  main  program  which  generate 
particular  outputs  when  flagged  to  do  so: 

TPDN  for  generating  output  of  the  height  profiles  of  temperature 
and  log^Q  of  pressure,  density  and  total  number  density  at 
given  height  increments  and  their  third  differences  with 
respect  to  height  (to  enable  the  smoothness  of  these 
profiles  to  be  examined  at  70  and  130  km  where  continuity 
in  the  second  height-derivative  has  been  formulated). 

CHPARM  for  generating  output  of  MSIS-86  diffusion  and  mixing  number 
densities  at  130  km,  i.e.  n^^.M^)  and  MSIS-86 

density  parameter  and  the  amended  value  (equation 
A5.15);  and  turbopause  height  z^  (equation  A. 96). 

CHNUMD  for  generating  output  of  log-^Q  of  the  number  density  of 
gas  constituents,  n(z,MJ,  at  given  height  increments  and 
their  third  differences  with  respect  to  height  (to  enable 
the  smoothness  of  these  profiles  to  be  examined  at  70  and 
130  km  where  continuity  in  the  second  height-derivative 
has  beeen  formulated). 

ATMOS  for  generating  output  of  temperature  and  of  log^Q  of 
pressure,  density,  total  number  density  and  individual 
gas  number  densities. 

ATMOSN  for  generating  output  of  temperature,  number  density, 
total  number,  pressure  and  density. 


Subprograms;  TEMPS,  P0LYZ12S,  BCZ12S,  BCATZ2S,  BCATZ2S,  SCATZ1S, 
EARTHS,  MSISINS,  LINES,  GTS5S,  CF1880S,  CFZ1Z2S. 

Input  files:  TPW12Q,  COEFFDQ,  WADJ(N1>D,  PARAMD. 

Output  file:  TABLEQ. 

TVAL1KMP  (TVALP  amended  to  give  1  to  interval  table  Program) .  TPDN  and 
CHNUMD  are  deleted  from  TVALP  and  changes  are  introduced  in 
conjunction  with  two  subprograms  that  are  added: 

TPDTB1  for  generating  output  of  separate  temperature ,  pressure 
or  density  tables -on  each  page  of  output  with  the  same 
format  and  1  km  height  interval  as  the  corresponding  tables 
in  the  report  'A  Global  Reference  Atmosphere  From  18  to  80 
km'  (Reference  1).  Such  tables  appear  in  Appendix  F. 

POWER  for  use  with  TPDTB1  to  calculate  the  power  of  10  needed 
in  the  pressure  and  density  tables. 

Subprograms  and  Input  files:  as  TVAL  (but  PRM86D1  replaces  PARAMD) . 
Output  file:  TAB86Q1 . 

TVAL5KMP  (TVALP  amended  to  give  5  KM  interval  table  Program) .  TVAL1KMP 
has  subprogram  TPDTB1  replaced  by  TPDTB5  and  a  few  associated 
changes : 

TPDTB5  for  generating  output  of  temperature,  pressure  and  density 
tables  with  a  5  km  height  interval  for  each  of  four 
selected  months  such  that  these  12  tables  fit  into  one 
page  of  output  as  shown  in  Appendix  F. 

Subprograms  and  Input  files:  as  TVAL  (but  PRM86D5  replaces  PARAMD). 
Output  file:  TAB86Q5. 
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APPENDIX 


OOEFFDQ 


TPW12Q 


BCATZ1S 


BCATZZS 


E  MEMO  LIST  OF  SUBPROGRAMS  AND  INPUT  DATAFILES 


(Coefficients  from  part  D  of  an  earlier  Q  (standing  for  output)) 
is  a  datafile  containing  the  coefficients  c  for  the  lower 
model  (Reference  1,  pp.  108-  109). 

(Temperature ,  Pressure  Waves  _1  and  2  from  an  earlier  2)  is  a 
datafile  containing  the  tables  of  amplitudes  and  phases  of  Ref.  1, 
pp.  110-  121,  which  define  the  longitudinal  dependences  of 
temperature,  pressure  and  density  for  the  lower  model.  (Only 
the  temperature  and  pressure  dependences  are  utilized). 

(Boundary  Conditions  at  height  Subprogram)  is 
SUBROUTINE  BCATZl(ZZ,DAY,MN,GLAT,GLONG,G) 
which  is  used  with  TDIFAVP4  (Appendix  D)  and  interpolates  coeff¬ 
icients  c  of  the  18-80  km  region  to  the  given  day,  DAY,  of 

I  lo 

month,  MN,  and  calculates  the  required  conditions  at  latitude,  GIAT, 
and  height,  ZZ  (=  =  70  km)  transferring  these  values  back  to  the 

calling  program  as  G(1),...G(4)  =  p(z^,M^  ),  100  W  ^  and  the  first 
two  derivatives  of  100  W  with  respect  to  height  at  height  z^ 
according  to  the  relations  in  Appendix  A2.  Longitude  (GLONG) 
dependence  is  included  unless  ILONG1  is  set  equal  to  zero  in  the 
main  program. 

(Boundary  Conditions  at  height  z^  or  at  lower  height  z  Subprogram)  is 
SUBROUTINE  BCATZZ (ZZ, DAY, MN.GLAT, GLONG, G) 
which  is  used  with  TVALP,  TVAL1KMP,  TVAL5KMP  (Appendix  D)  and  is 
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BCATZ1S(ZZ,MY,MN,GIAT,GL0NG,G)  with  an  additional  section  of 

instructions  such  that  when  ZZ  =  70  it  provides  the  same  G(l),... 

G(4)  as  BCATZ1S  and  for  other  ZZ  (=  z)  it  provides  G(l)  = 

ptz,^  )  and  G(2)  =  100 W  *  at  height  z  (less  than  70  km). 
w2 

GTS5S  is  the  MSIS-86/CIRA  1986  neutral  atmosphere  model  of  15  March 

1986  prepared  by  A.E.  Hedin  as 

SUBROUTINE  GTS5(IYD,SEC,ALT,GIAT,GL0NG,STL,F107A,F107,AP, 

MASS,D,T) 

with  the  additional  facility  to  transfer  parameters  to  other 

programs  through  the  common  block  AA  defined  by 

COMMON/ AA/  DDF(8) ,DMX(8) ,HC04,HC16,BLKl,HC32,HC40,BLK2,HC01, 
HC14 , ZC04 , ZC16 , BLK3 , ZC32 , ZC40 , BLK4 , ZC01 , ZC14 , BLK5 , 
RC16 ,BLK6 ,BLK7 ,BLK8 ,BLK9 ,RC01 ,RC14,BLK10,HCC16 , 

BLK1 1 ,  BLK1 2 ,  BLK1 3 ,  BLK14 ,  HCC01  ,HCC14 ,  BLK1 5 ,  ZCC1 6 , 
BLK16 ,BLK17 ,BLK18 , BLK19 , ZCC01 ,ZCC14 ,RCU( 8 ) 

where  DDF(I),  DMX(I)  are  respectively  the  diffusion  profile  number 

density  nj(z,Mj)  and  the  mixing  profile  number  density  n^z.Mj) 

corresponding  to  I  =  1,2, 3,4, 5, 7, 8.  RCU(I)  are  likewise  the 

values  of  R^  (Appendix  A5). 

BCATZ2S  (Boundary  Conditions  at  height  Subprogram)  is 

SUBROUTINE  BCATZ2(IYD,SEC,Z2,GLAT,GL0NG,STL,F107A,F107,AP,G) 

which  calls  GTS5  to  calculate  the  MSIS-86  conditions  required  at 

height  z2  (=  130  km)  and  transfer  them  to  the  calling  program  as 

G(1),...G(6)  =  Pm(z2,MN  Pd^z2’^N  ^  and  t^ie  first  three 

_1  2  2 

derivatives  of  100  W  with  respect  to  height  at  height  z2  accord¬ 
ing  to  the  relations  in  Appendix  A3.  (G(6)  is  not  utilized). 

BCZ12S  (Boundary  Conditions  at  z^  and  z2  Subprogram)  is 

SUBROUTINE  BCZ12 ( Z1 , Z2 , IDAY ,MN ,GLAT,GLONG , SLT , F107A ,F107 , 

AP ,G1 ,G2 ) 


and  is  the  calling  program  of  (1)  either  BCATZl  (if  used  with 
TDIFAVP4)  or  BCATZZ  (if  used  with  TVALP,  TVAUKMP  or  TVAL5KMP) 
to  obtain  G(I)  as  G1(I)  to  which  the  sunspot  number  dependent 
changes  are  applied  by  calling  SCATZ1  and  (2)  BCATZ2  to  obtain 
G(I)  as  G2(I). 

SCATZ1S  (Solar  Cycle  at  height  Subprogram)  is 

SUBROUTINE  SCATZ1(ZZ,IYD,GLAT,RNDEL,GDEL) 
which  calculates  incremental  values  GDEL(  I )  arising  from  an 
incremental  change  of  sunspot  number  RNDEL  according  to  the 
relations  (A2.22)  and  (A4.15). 

RZDATA  (Sunspot  number,  R^,  Datafile)  contains  monthly  mean  sunspot  numbers 
from  January  1957  to  January  1972,  the  interval  of  time  within  which 
the  launch  dates  of  single  rocket  profiles  fall. 

P0LYZ12S  (Polynomial  coefficients  for  height  range  z-^  to  Subprogram)  is 

SUBROUTINE  P0LYZ12 ( Z1 , Z2 , IDAY ,MN,GLAT ,GL0NG ,SLT,F107A,F107 , AP,B) 
which  calculates  the  polynomial  coefficients  b  from  (A1.15)  having 
first  obtained  £_  which  involves  the  relations  of  Appendix  A4  to 
obtain  from  (A4.17). 

EARTHS  (Earth  radius  and  gravity  Subprogram)  is 
SUBROUTINE  EARTH (ALAT,GEPHI  ,RPHI) 
which  calculates  g^,  and  r ^  from  (A. 8)  and  (A. 9)  at  latitude 
ALAT  (*f). 

MSISINS  (MSIS  variations  included  Subprogram)  is 
SUBROUTINE  MSISIN 

which  lists  those  of  the  23  variations  of  MSIS-86  that  are  being 


omitted  when  the  Subroutine  GTS5  is  called. 

THADATA  (Temperature  high  altitude  Datafile)  contains  observed  temperatures 
at  5  km  height  intervals  from  70  to  130  km  at  different  sites. 

POLYTDD  (Specifies  polynomial  for  obtaining  temperature  differences  Datafile) 
is  a  file  of  input  parameters  for  TDIFAVP4:  IH1,  IH2  =  the  range 
of  data  in  THADATA  to  be  utilized,  NOBS  =  the  number  of  single 
profiles  to  which  a  monthly  mean  is  equated  for  purposes  of 
weighting  (finally  taken  as  3);  ISLT  =  0  or  1  according  to  whether 
local  solar  time  is  excluded  or  included  in  MSIS-86,  IL0NG1  = 

ILONG2  =  0  or  1  according  to  whether  longitude  dependences  are 
excluded  or  included  in  lower  model  or  upper  model  (i.e.  MSIS-86), 
NCOND  =  number  of  coefficients  (=  7)  of  the  interpolating  polynomial 
between  heights  and  NF107  (=  120)  and  NP  (=  10)  are  values 
of  solar  activity  parameters  F^q  j  and  which  are  adopted  when 
no  other  values  are  specified,  NRUN  =  reference  number  assigned  to 
a  particular  computer  run,  IMSIS  =  0  for  normal  use  of  TDIFAVP4 
and  =  99  to  give  tables  of  T  -  TMSIS  where  T  is  a  polynomial  model 
temperature  and  IMSIS  is  the  corresponding  MSIS-86  value,  ISAX 
(=  7)  and  INAX  (=  2)  are  the  values  of  S  and  N  in  Appendix  B, 

NOSCZ1  *  0  for  RNDEL  to  be  set  zero  in  BCZ12  so  that  no  sunspot 
number  dependence  is  introduced  into  the  lower  model. 

DtMWAXD  (Dunmy  set  of  a  coefficients  Datafile)  is  datafile  of  a  all 
■  — J  — sn  —  sn 

of  which  are  zero. 

TEMPS  (Temperature  etc  Subprogram)  is 

SUBROUTINE  TEMP(ZHT1,ZHT2,ZHTD,IHEND,Z1,Z2,IDAY,MN,GLAT,GLJ0NG, 

SLT ,  F107A ,  F107 ,  AP  ,TT ,  DN ,  PR ,  TOTND  ,RCL ,  ZZTURB ,  DDFZ2 ,  DHXZ2 ) 


LINES 


CFI880S 


CFZ1Z2S 


PARAMD 


which  calculates  at  each  of  IHEND  heights  from  ZHT1  to  ZHT2  at  a 
height  interval  of  ZHTD  the  temperature  TT,  the  number  densities 
of  individual  gas  constituents  DN(I),  1^  6.  and  the  total  mass 
density  DN(6),  the  total  no.  density  TOTND,  the  parameter  RCL(I) 

(=  from  equation  A5.16),  the  ' turbopause '  heights  ZZTURB(I) 

(=  of  equation  A. 97),  diffusion  profile  number  densities 
DDFZ2(I)  (=  n^(z2#Mj))  and  mixing  profile  number  densities  DMXZ2(I) 
(=  nm(z2,Mj))  at  height  Z2  of  the  gas  constituents. 

(Line  Subprogram)  is 

SUBROUTINE  LINE (CHAR) 

which  reads  in  and  writes  out  a  line  of  characters. 

(Coefficients  for  18-  80  km  model  Subprogram)  is 
SUBROUTINE  CF1 880 (CC, DELS, WAV) 

which  reads  in  TPW12Q  and  derives  DELG  and  WAV  for  use  in  BCATZZS 
and  then  reads  in  COEFFDQ  into  CC. 

(Coefficients  agn  for  z^  to  model  Subprogram)  is 
SUBROUTINE  CFZ1Z2 

which  reads  in  and  writes  out  a  . 

sn 

(Parameter  Datafile)  runs  TVALP  with  ZHT1,  ZHT2,  ZHTD; 

MN1,  MN2,  MND;  LAT1,  LAT2,  LATD;  LNG1,  LNG2,  LNGD;  ISLTl, 

ISLT2 ,  ISLTD  as  the  ranges  of  height,  month,  latitude,  longitude  and 
local  solar  time  and  their  respective  increments  at  which  values 
are  required  to  be  evaluated;  IXLONG  =  0  gives  zonal  mean  values; 
ICHPAR,  ITPDN,  ICHNUM,  IATMOS,  ITMOSN  =  0,  1  for  output 
subroutines  CHPARM,  TPDN,  CHNUM,  ATMOS  and  ATMOSN  not  to  be 
called  or  to  be  called;  IXSLT  =  0  gives  diurnal  mean  values 
in  MSIS-86;  solar  activity  parameters  for  which  TVALP 
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PRM86D1 


PRM86D5 


is  to  run  are  listed  as  F107A,  F107,  AP(1) , .  .AP(7) ,  SW(9),  and 
N0SCZ1,  where  SW(9)  controls  the  use  of  AP(I)  in  GTS5  (i.e. 

MSIS-86)  and  N0SCZ1  =  0, 1  excludes  or  includes  solar  cycle 
dependence  in  the  lower  model  according  to  SCATZ1. 

(Parameters  for  MS  IS -86  Datafile  with  TVAL1KMP)  gives  ZHT1  (=65.), 

ZHT2  (=  135. ) ,  ZHTD  (=  1. ) ;  .  LATl(=-80);  LAT2(=80);  LAID  (=10); 

LNG1,  LNG2,  LNGD;  ISLT1,  ISLT2,  ISLTD  as  the  ranges  of  height, 
latitude,  longitude  and  local  solar  time  and  their  respective 
increments  at  which  values  are  required  to  be  evaluated; 

IXLONG  =  0  gives  zonal  mean  values;  IXSLT  =  0  gives  diurnal  mean 
values  in  MSIS-86,  ICHPAR(=0),  IATMOS(=0);  ITM0SN(=0), 

ITPDTB  (=  1)  .  for  output  subroutines  CHPARM,  ATMOS,  ATMOSN,  TPDTB1 
not  to  be  called  or  to  be  called;  IMSIS  =  0  for  normal  use,  other¬ 
wise  tables  of  values  of  MSIS-86  are  generated;  tables  are 
generated  for  given  dates  specified  by  day  of  month  and  month  as 
IDAT(I),  MON(I),  I  =  1,..NDAYS;  IDAT(I) =  0  gives  mid-month  values; 

IXMN=0  gives  annual,  mean,  =1  for  given  dates  as  specified  above; 
solar  activity  parameters  for  which  TVAL1KMP  is  run  and  their 
controlling  parameters  SW(9)  and  N0SCZ1  are  supplied  as  in  PARAMD. 
(Parameters  for  MSIS-86  Datafile  with  TVAL5KMP)  is  the  same  as 
PRM86D1  with  different  height  and  latitude  ranges;  ZHT1  (=70.), 

ZHT2  (=  130.),  ZHTD  (=  5.);  LATl(=-80),  LAT2(=80),  LATD(=20). 

< 


i 
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APPENDIX  F 


EXAMPLES  OF  TABULATIONS  OF  TEMPERATURES , 
PRESSURES  AND  DENSITIES 


Two  formats  have  been  devised  for  tabulating  diurnal  and  zonal  means  of 
mid-month  values  of  temperature,  pressure  and  density: 

(i)  With  a  1  km  height  interval  and  a  10°  increment  of  latitude 
such  that  each  of  the  three  parameters  requires  1  page  per 
month.  The  height  range  is  65  - 135  km,  the  values  for 
65-  70  km  being  those  of  the  lower  model  and  those  for 
130-  135  km  being  from  the  upper  model. 

(ii)  With  a  5  km  height  interval  and  a  20°  increment  of  latitude 
such  that  all  three  parameters  for  4  consecutive  months  fit 
on  one  page  for  the  height  range  70  -  130  km.  The  months 
on  each  page  are  grouped  as  January  -  April,  May- August, 
September  -  December . 

The  examples  presented  are  all  for  solar  activity  parameters  F^g  ^=150, 

A  =4  and  for  Case  2  of  Section  4.2. 

P 
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